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1. Probing dense matter of elementary particles 



Ultra-relativistic heavy ion collisions offer the unique opportunity to probe highly 
excited dense nuclear matter under controlled laboratory conditions. The compelling 
driving force for such studies is the expectation that an entirely new form of matter 
may be created from such reactions. That form of matter, called the Quark Gluon 
Plasma (QGP), is the QCD analogue of the plasma phase of ordinary atomic matter. 
However, unlike such ordinary plasmas, the deconfined quanta of a QGP are not directly 
observable because of the fundamental confining property of the physical QCD vacuum. 
What is observable are hadronic and leptonic residues of the transient QGP state. 
There is a large variety of such individual probes. Leptonic probes, 7, e + e~, p^ pr 
carry information about the spectrum of electromagnetic current fluctuations in the 
QGP state; the abundance of quarkonia T, T' (also observed via carry 

information about the chromoelectric field fluctuations in the QGP. The arsenal of 
hadronic probes, tt, K, p, p, A, S, Q, 0, p, ... provide information on the quark flavor 
chemistry and baryon number transport. Theory suggests that with decays such as 
p — ► e + e~ the properties of the hadronization and chiral symmetry breaking can 
be indirectly studied. Quantum statistical interference patterns in mr, KK, pp, AA 
correlations provide somewhat cloudy lenses with which the space-time geometry of 
hadronic ashes of the QGP can be viewed. The detailed rapidity and transverse 
momentum spectra of hadrons provide barometric information of pressure gradients 
during the explosive expansion of the QGP drop. 

The central problem with all the above probes is precisely that they are all indirect 
messengers. If we could see free quarks and gluons (as in ordinary plasmas) it would be 
trivial to verify the QCD prediction of the QGP state. However, nature choses to hide 
those constituents within the confines of color neutral composite many body systems - 
hadrons. 

The QGP state formed in nuclear collisions is a transient rearrangement of the 
correlations among quarks and gluons contained in the incident baryons into a larger 
but globally still color neutral system with however remarkable theoretical properties. 
The task with heavy ion reactions is to provide experimental information on that 
fundamental prediction of the Standard Model. 

This topical review covers current (1998) theoretical and experimental attempts 
to disentangle popular scenarios on QGP signatures from the complex, off-equilibrium 
physics. The start of the RHIC experiment program is only a year away. This will be 
a dedicated machine to study the QGP. Nevertheless a very large effort has been made 
at the AGS and SPS over the last 12 years and has resulted in an impressive amount 
of exciting new findings. The search for the QGP can be traced via the proceedings of 
the High Energy Heavy Ion Studies § g g § g g 0, g g and Schools |TT], |TJ, |HJ 
and the "Quark Matter", "Nucleus- Nucleus" and "Strange Quark Matter" conferences 

n m m 0, m m m n n, m m n, m n m he m m h m ■ s ° me textbooks 



34] , j35| , 361, WA EH) an d a vast number of review articles have been published - reference 
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samples of early |9|, fl(J 0, |g §g (45], || and of the latest g7|, f4|, g|, 0, || review 



papers are given here. The list of more than 500 references given in this topical review 
is by no means complete. Apologies are offered to those whose contributions could not 
be included into the references. 

2. QCD matter and relativistic heavy-ion collisions 

2.1. Infinite stationary systems in equilibrium 

The deconfinement phase transition and chiral symmetry restoration 
Phase transitions are among the most dramatic many body effects in physics. Examples 
for restored symmetry via a phase transition at high temperatures, Tc, are ferro- 
magnetism, super-conductivity and the solid - liquid phase transition. In nuclear 
physics evidence for a liquid gas phase transition of nuclear matter has been claimed for 
temperatures of T pa 5 MeV ||52j| . Phase transitions to abnormal nuclear matter states 



at high densities have also been predicted early on |53], p4| . 

QCD is a non-abelian gauge theory, it's basic constituents are quarks and anti- 
quarks interacting through the exchange of color-charged gluons. At short space- 
time-intervals - large momentum transfers - the effective coupling constant decreases 
logarithmically ( "asymptotic freedom" - meaning weak coupling of quarks and gluons) 
whereas it becomes strong for large distances and small relative momenta. This results 
in the phenomena of chiral symmetry breaking and quark- gluon confinement. 

At very high temperatures and densities, in the domain of weak coupling between 
quarks and gluons, long range interactions are dynamically screened |)5[ 56fl . Quarks 



and gluons are then no longer confined to bound hadronic states ("deconfinement"). 
Furthermore, chiral symmetry is restored - for baryon-free matter - apparently at the 
same temperature T c . This novel phase of nuclear matter is called the quark gluon 



plasma 55 



A transition from the deconfined quark-gluon phase to confined color singlet states 
has (probably) occured during the rapid expansion of the early universe. Temperatures 
were very high then, but the net baryon density was small. Therefore one often assumes 
zero baryon chemical potentials in calculating the thermodynamic properties of strongly 
interacting matter in the early universe. It is sought to re-establish these conditions and 
thus enable a study of quark deconfinement in the laboratory via heavy ion collisions 
|57| , |58fl . At the highest in the near future obtainable energies (LHC) the initial net 
baryon density may be around nuclear matter ground state density. The entropy 
per baryon ratio, however, is estimated to be in the order of 10 3 to 10 4 |)]| (early 
universe: 10 9 ) and thus a vanishing baryon chemical potential is considered a viable 
approximation. 

The energy densities currently achievable in ultra-relativistic heavy ion collisions 



at the AGS and SPS are on the order of 0.5-10 GeV/fm |60| and temperatures are in 



the range of 100-200 MeV. At temperatures Tq ~ 150 — 200 MeV the effective coupling 
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constant of QCD, however, is still on the order of 1. Therefore, perturbative techniques 
of QCD [62|, |63| are not applicable. The situation may change at extremely high 
energies (LHC?) where QCD predicts a large cross section for minijet production. Then 
there may exist a possibility to reach a new regime of large parton densities at a small 



coupling constant. The QCD interaction in this regime would be highly nonlinear |64 
Recent calculations of the Debye screening mass using lattice gauge simulations, however, 
indicate that perturbative QCD techniques may not even be applicable at such energies 



Perturbation theory 

The fundamental asymptotic property of QCD leads to the naive expectation that 
the properties of a QGP can be calculated via perturbation theory. However, due to 
infrared divergences, especially in the chromo-magnetic sector, perturbation theory may 
not even be applicable for temperatures far above Tq. QCD perturbation theory has 
been improved |)2j by resummation to screen color-electric divergences. These "hard- 
thermal-loop" methods, originally developed for zero baryon density, have now been 
extended to finite densities [|BEf| . Nevertheless, its applicability at temperatures and 
densities accessible to experiment is severely limited. At energy densities on the order 
of 20 GeV/fm 3 the temperature is in the order of 600 MeV and the coupling constant 
g is still of the order of 2, which invalidates the necessary assumption gT <C T for the 
hard-thermal-loop resummation scheme. Only at the Planck scale is QCD a weakly 
interacting theory as QED |5l|. Recent calculations of the perturbative contributions 



up to 0(g 5 ) |)7j to the pressure are shown in fig.|T[ The oscillation of the results suggest 



a zero radius of convergence of thermal pQCD. 
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Figure 1. Perturbative contributions up to 0(g 5 ) J67J to the pressure vs. the 
coupling constant. 
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Lattice Gauge Theory 

Lattice gauge simulations of QCD ^{J, provide therefore the only rigorous method 
to compute the equation of state of strongly interacting elementary particle matter. In 
principle both, the non-perturbative hadronic matter and the non-perturbative QGP 
phases of QCD can be investigated. The main disadvantage of lattice simulations is 
the practical restriction to finite, periodic, baryon free systems in global equilibrium, a 
scenario far from the highly inhomogeneous off-equilibrium situation found in complex 
heavy-ion reactions. Technically, the strong dependence of the results on the lattice 
spacing and periodic box size is presently a problem. Nevertheless, lattice data provide 
the most compelling theoretical evidence of a rapid transition region from the confined 
to the QGP state. 

Lattice calculations allow at least for the computation of thermodynamic averages of 
different quantities related to hadron masses and to the phase transition (in the infinite 
volume and zero baryon number limit). There have been considerable improvements in 



algorithms |70[ and in computing power in recent years. For finite temperature full QCD 
simulations lattices with spatial sizes of 24 3 and 48 points in Euclidean time direction 
have been used |7T|, [72|], while for pure gauge theory (without quarks) lattices of 32 3 x 12 

m 
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Figure 2. Left: Lattice calculation (circles) of the equation of state for two flavor 
QCD |7^| . The lines show fits using 0(4) scaling and extrapolations to zero quark 
mass. The critical temperature is in the order of 140 MeV. Right: schematic overview 
of theoretical phases of QCD matter. 



Lattice calculations yield a critical temperature of Tq = 265^5° MeV in the 
quenched approximation |7l| - where neither dynamical quarks, nor a chiral phase 
transition exist. Simulations including dynamical quarks at \lb — indicate a critical 
temperature in the order of Tc = 140 MeV (see figure |2|). However, in this case finite 
size effects of the lattice have not yet been fully overcome and the precision is not 
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as high as in the quenched case J7[J. The inclusion of the second most important 
thermodynamic variable, the chemical potential /i# into a full fletched 1QCD calculation 
is, presently, still out of reach. This raises a practical question, whether conclusions 
based on \xb = estimates, might misguide physical argumentation for observables in 
nuclear collisions. This warning is particularly appropriate for those QGP-signals, where 
a 50% quantitative change of an observable is used to differentiate QGP production 
scenarios from ordinary hadronic transport ones. 

The behavior of order parameters as a function of temperature, such as the quark 
condensate, indicate that the transition between the confined and deconfined state 
of QCD may show a discontinuity of some thermodynamic derivative - i.e. a phase 
transition. The order of this phase transition is crucial for some proposed signatures 
of the QGP. Many striking signatures depend heavily on the assumption of a first 
order phase transition and the existence of a mixed phase of QCD matter. For pure 
SU(3) gauge theory and for full QCD with four massless flavors of dynamical quarks, 
lattice QCD results and universality class arguments |76|, [77| indicate a first order phase 
transition. For two (massless) flavors, however, these arguments predict ]F7], [78| that if 
the phase transition is of second order, it should have the critical exponents of the 0(4) 
Heisenberg anti-ferromagnetic spin-model in three dimensions. Numerical evidence for 
this was recently obtained J79|, suggesting that the transition in this case is indeed 
of second order. With three light flavors the chiral phase transition corresponds to a 
change in a continuous symmetry (i.e. chiral symmetry) and is automatically of second 
order. 

For the most realistic case of QCD with two flavors of light quarks with masses 
between 5 and 10 MeV and one flavor with a mass around 200 MeV, the situation 
remains unclear: the order of the phase transition seems to depend on the numerical 
values for the masses of the light and heavy quarks [80]. If the latter is too heavy, the 
transition might be smeared out to a mere rapid increase of the energy density over 
a small temperature interval. In this case the use of simple deconfinement scenarios 
may lead to wrong expectations for observables. The elementary excitations in such a 
phase transition region ought not be described by quarks and gluons but could physically 
resemble more hadronic excitations with strongly modified "in-medium" properties |8l[ . 

In any case the QCD phase transition observed on the lattice is - when dynamical 
quarks are included - only very weakly first order, if there is a discontinuity at all. The 
latent heat across the discontinuity is at most a small fraction of the total jump in the 
normalized entropy density s/T 3 between the hadronic phase and the asymptotic QGP. 
A real phase-coexistence region between hadrons and (possibly strongly interacting, 
but deconfined) quarks and gluons, as often discussed in the past, seems therefore no 
longer a realistic possibility. However, nothing is known from first principles about the 
phase transition at finite baryon density and therefore a strong first order transition 
(with a phase coexistence region) is still a possibility at large baryon densities [^, Q . 
An additional complication is that for systems of finite volume (V< 100 fm 3 ) the 
deconfinement cannot be complete. Fluctuations lead to a finite probability of the 
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hadronic phase above T C - The sharp discontinuity (e.g. e/T A ) is thus smeared out |84 

Purely hadronic models, such as the a — uj Model or the linear a-model exhibit 
a similar phase transition (from normal to abnormal nuclear matter), but are not 
constrained to ps = 0. The equation of state for nuclear matter does not only depend 
on temperature and density but will also depend on the net-strangeness content, which 
may be non-zero in subsystems (e.g. individual phases) present in heavy ion collisions. 
A schematic view of the resulting complex multidimensional phase-diagram is illustrated 
on the right in figure ^. 



2.2. N on- equilibrium models 

In order to connect the theoretical thermodynamic properties of a QGP with 
experimental data on finite nuclear collisions, many non-equilibrium dynamical effects 
must also be estimated. Transport theory is the basic tool to address such problems. 
Non-equilibrium effects are certain to arise from the rapid time-dependence of the 
system (even the use of the term "state" seems questionable), finite size effects, 
inhomogeneity, iV-body phase space, particle/resonance production and freeze-out and 
collective dynamics. Such microscopic and macroscopic (hydrodynamical) models 
attempt to describe the full time-evolution from an assumed initial state of the heavy 
ion reaction (i.e. the two colliding nuclei) up to the freeze-out of all initial and produced 
particles after the reaction. Hydrodynamical models neglect most of these effects 
by making the assumption that the initial condition can be assumed to be in local 
thermal equilibrium and that local equilibrium is maintained during evolution. Fireball 
models simply parameterize final spectra and abundances via freeze-out parameters, e.g. 
T^BjVf. However, the initial condition in nuclear collisions is a coherent state \AB) 
of two quantal (T = 0) nuclear systems. A non-equilibrium quantum evolution of \AB) 
introduces complex high order Fock-State components. A key dynamical assumption 
is that decoherence occurs rapidly during the early phase of the collision yielding a 
mixed state density matrix (with S = Trplnp > 0). There is no theorem to insure that 
p evolves to a local equilibrium form exp(— w M p M /T) at any time during the reaction. 
That can only be tested via a transport theory approximation to the evolution equations. 
The question of the form of the initial state p(t ) must still be addressed, but once 
that is specified, transport theory can reveal if local equilibrium is achieved and what 
observables are least sensitive to uncertainties in p(tq). 

Depending on the most convenient basis for expanding p(r ), transport theory 
assumes different forms. At low energies the initial ensemble is most conveniently 
described in terms of mesons and baryons. Here hadronic transport theory is 
appropriate. At collider energies, pQCD minijet processes are expected to produce 
a high density mostly gluonic gas. In that regime parton cascade models are more 
appropriate. 
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Parton cascades 
Parton cascade models 



5], [87], ^SJ, ^9|, |9(| evolve partonic degrees of freedom. They 
are therefore mostly applied to study the initial compressional and the high density 
phase of ultra- relativistic heavy ion collisions (collider energies, y/s > 200 GeV). These 
models all contain the general structure p6| : 

(i) Initialization: the nucleons of the colliding nuclei are resolved into their parton 
substructure according to the measured nucleon structure functions and yield the 
initial parton distributions. 

(ii) Interaction: parton interactions as described by perturbative QCD are used to 
model the evolution of the ensemble of partons during the course of the collision. 
This includes multiple scatterings together with associated space-like and time- 
like parton emission processes before and after each scattering. The sequence of 
scatterings is, however, incoherent and the neglect of quantum interference effects 
is questionable. 

(iii) Hadronization: partons are recombined or converted via string fragmentation into 
final hadron states. 

The propagation is performed on straight lines - soft non-perturbative collective field 
effects have so far been neglected. On the other hand, hadronization has to be modeled 
by brute force to mock up confinement in the final reaction stage. 
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Figure 3. Comparison of time evolution of transverse energy per rapidity in 
analytic kinetic theory results with numerical parton cascade calculations ]92| ] . Strong 
deviations from hydrodynamic behavior are visible. 



One of the central issues addressed by parton cascades is the question of energy 
deposition processes in space-time as well as momentum space. Partonic cascades 
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predict that roughly 50% of the expected energy deposition at RHIC and a larger 



fraction at LHC takes place at the partonic level [91]. Rapid thermalization is caused 



by radiative energy degradation and spatial separation of partons with widely different 
rapidities due to free streaming; transverse momentum distributions of initially scattered 
partons are almost exponential if radiative corrections are taken into account ||49|| . For 
RHIC energies thermalization is predicted on a proper time scale of 0.3 - 0.5 fm/c ||91|| . 
A recent analysis of parton cascade evolution |)2| shows that local equilibrium is not 
maintained due to rapid expansion. Very large dissipative corrections to hydrodynamics 
appear. The thermalized QGP is initially gluon rich and depleted of quarks due to the 
larger cross section and higher branching ratios for gluons ||93|1 . Chemical equilibrium 



is achieved over a time of several fm/c [94), f)5]. This may be reduced if higher order 



pQCD processes are taken into account |96] 



Hadronic transport models 

Hadronic transport models treat relativistic heavy-ion collisions as sequences of 
binary/ iV-body collisions of mesons, baryons, strings and their constituents, diquarks 
and quarks. The real part of the interaction can be obtained in principle from G-Matrix 
calculations, with the in-medium self-energy and the imaginary part is modeled via hard 
scattering cross sections |97|, |9|, |9|, |10g, [101], [102], [103], [104|, [TOol |T06l, g07[ TO [109|, [TT0[ . 
For high beam energies most models include particle production via string formation 
- either using the Lund [ 1 IT , 112 , 113 | or a pomeron exchange scheme | 114 |. Partonic 
degrees of freedom are not treated explicitly and therefore these models do not include 
a phase transition. However, some models contain further speculative scenarios such as 
color-ropes |115| , |116|| , breaking of multiple-strings [|11711 or decay of multi-quark droplets 



118|| which clearly go beyond hadronic physics. 

Hadronic transport models are critical for assessing the influence of ordinary or 
exotic hadronic phenomena on the observables proposed to search for a QGP. They 
therefore provide a background basis to evaluate whether an observable shows evidence 
for non-hadron physics. 



Nuclear fluid dynamics 

NFD is so far the only dynamical model in which a phase transition can explicitly be 
incorporated (see e.g. pD|, |119j , |120|| for details). This is possible since the equation of 
state (including a phase transition) is a direct input for the calculations. However, NFD 
is an idealized continuum description based on local equilibrium and energy-momentum 
conservation. Therefore it is very well suited to study kinematic observables such as 
collective flow. Since NFD is a macroscopic kinetic theory it is not directly applicable 
to the study of hadron abundances and particle production. However, NFD calculations 
predict (local) temperatures and chemical potentials which can be used, e.g. by chemical 
equilibrium calculations of hadron abundances, to study particle production. Different 
observables predicted by nuclear fluid dynamics will be discussed in section |3~3| . 

In the ideal fluid approximation (i.e. neglecting off-equilibrium effects), the EoS is 
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the only input to the equations of motion that relates directly to properties of the matter 
under consideration. The EoS influences the dynamical evolution of the system, and final 
results are uniquely determined. The initial condition can be chosen from two colliding 
nuclei (in a full 3D calculation with up to three fluids) or an equilibrated QGP or 
hadronic matter with prescribed temperature and chemical potential and velocity/flow 
profiles (for simpler, more schematic calculations). The time-evolution is then studied 
until hadronic freeze-out for which a decoupling (freeze-out) hyper-surface needs to be 
specified. 

However, the ideal fluid ansatz is only a rough approximation. In the parton cascade 
study [|92] for example, large deviations from even the Navier Stokes fluid approach were 
found. 



3. Observables: prospects and limitations 

As we have seen in the previous section it is obviously difficult to find a robust theoretical 
description of relativistic heavy ion collisions involving the QCD phase-transition to 
predict observables. Not only is even the order of the phase-transition from fi = 
not known from the ab-initio 1QCD calculations, but also has the physical situation 
of present or near future relativistic heavy ion collisions, namely finite not been 
addressed yet in this theory. However, even if this would be the case, one would 
only know the behavior for static infinite systems. The second major unknown is 
the influence of the non-equilibrium evolution on the (small) many-body system. The 
very nature of even the thermodynamic limit of a QGP is not completely understood. 
Real time response has only been studied via pQCD, which however may have zero 
radius of convergence in g in the thermodynamic limit. Theory in this situation 
can thus serve mainly to motivate particular experimental studies and provide overall 
consistency checks in the interpretation of data. Data are needed to fix the uncertain 
phenomenological parameters of the transport models, while such model calculations 
with plausible parameters are essential to motivate the taking of the data in the first 
place. This symbiotic relation between theory and experiment in this field is very 



important as emphasized also for example by Van Hove ||121|| and Kajantie ||122 . 

A strategy for the detection of quark matter - in our opinion - must collect at 
least circumstantial evidence from several "signals" or anomalies. In the following we 
discuss each of the individual signals. The strategy does then consist in a systematic 
variation of an external parameter (system size, impact parameter and - in particular 
- bombarding energy); i.e. the measurement of excitation functions of several signals 
which in the case of a phase-transition show simultaneously the predicted anomalous 
behavior. 



11 



3.1. Creation of high baryon density matter: nuclear stopping power 
Theoretical concepts 

It has been proposed more than two decades ago that head-on collision of two nuclei can 



be used to create highly excited nuclear matter |57|, |58| - The longitudinal momentum 
is converted via multiple collisions into transverse momentum and secondary particles, 
causing the creation of a zone of high energy density. Nuclear shock waves have been 
suggested as a primary mechanism of creating high energy densities in collisions with 



^fs < 20 GeV. |57], [58], 54]]. This is analogous to the well known Rankine-Hugoniot 
analysis of ordinary dense matter up to ~ 1 Mbar pressures. In the nuclear shock 
wave case, the Rankine-Hugoniot analysis predicts that pressures up to 10 23 Mbar 
~ 100 MeV/fm 3 may be reached. 

The term nuclear stopping power ||123|| characterizes the degree of stopping which 



an incident nucleon suffers when it collides with another nucleus. For A+A collisions 
stopping manifests itself in a shift of the rapidity-distributions of the incident nucleons 
towards mid-rapidity. The heaviest systems available, such as Pb+Pb or Au+Au, are 
best suited for the creation of high baryon density matter. 

The shape of the baryon rapidity distribution can give clear indications on the onset 
of critical phenomena: Due to the strong dependence of the baryon rapidity distribution 
on the baryon-baryon cross section ||124| , |125| , |126|| , a rapid change in the shape of the 



scaled dN/d(y/y p ) distribution with varying incident beam energy is a clear signal for 
new degrees of freedom which show up during the reaction (i.e. deconfmement), e.g. due 



to phenomena such as critical scattering ||127|| . The width of the dN / d(y / y p ) distribution 



for baryons is inversely proportional to their cross section. 

Hadronic transport model calculations have predicted stopping for heavy collision 
systems at CERN/SPS energies ||128| , |129|| (see figure |]). Even for RHIC energies the 



central rapidity zone is not expected to be net-baryon free. RQMD has predicted a 
net-baryonnumber density of > 10 at mid-rapidity [ |130|| and HIJING/B yields similar 
estimates [|131|| (see figure^). 



The creation of a zone of high baryon and energy density around mid-rapidity 
results in massive excitation of the incident nucleons. A state of high density resonance 



matter may be formed |58|, |135| , |136|| . Transport model calculations indicate that this 
excited state of baryonic matter is dominated by the A 12 32 resonance. They predict a 
long apparent lifetime (> 10 fm/c) and a rather large volume (several hundred fm 3 ) for 
this A— matter state in central Au+Au collisions at the AGS ||137|| (see figure 0). 

The degree of stopping can furthermore be used to estimate the achieved 
energy density in the course of the collision within the Bjorken scenario of scaling 
hydrodynamics For such an estimate not the rapidity distribution of the incident, 
leading particles is required, but that of secondary particles, those produced during the 
course of the reaction. One often assumes that particles produced at y = ycu originate 
from the central reaction zone at z = and the initial proper time r . The rapidity 
distribution of these produced particles could then be used to estimate the initial energy 
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Figure 4. Left: RQMD prediction [ 132 1 of stopping in central Au+Au collisions at 
10.6 GeV/nucleon. The preliminary data are from the E866 collaboration [133]. Note 
that the current status of data analysis indicates a flatter shape for the experimental 
distribution [ 134 1 . Right: RQMD prediction |12£] of stopping in central Pb+Pb 
collisions at 160 GeV/nucleon. 



density in the central reaction zone: 
rriT dN 
u a y=ycM 

Here A is the transverse overlapping region area in the collision and m-r the transverse 
mass of the produced particles. The proper production time r is very uncertain and 
estimates are on the order of 0.5 - 1 fm/c. Estimates for the CERN/SPS energy density 
at proper time To ~ 1 fm/c are on the order of e ~ 1 to 5 GeV/fm 3 ||60|| , with baryon 
densities up to p < 1 fm~ 3 . In ref. [|13l|| extrapolations to RHIC suggest that energy 
densities up to 20 GeV/fm 3 at p ~ 2p may be reached (see figure |5|). 



Experimental status 

At AGS and SPS an extensive investigation of the nuclear stopping power is near 
completion. Proton-proton [ |139| | and peripheral nucleus-nucleus interactions at AGS 
|140| , |141|| and SPS [|142|1 energies yield a forward-backward peaked dN/dy distribution 
in the CM. frame, and a low degree of baryon stopping. 

A higher degree of stopping is observed for central collisions of intermediate 
mass nuclei (Si+Si at AGS, S+S at SPS): The rapidity distribution is flat at CM. 
rapidities, two broad bumps are observed between projectile/target and CM. rapidities 
respectively ||140| , |141| , |142|| . The heaviest collision systems (gold and lead respectively) 
exhibit the largest stopping power and thus correspond to the creation of the highest 
baryon densities: At AGS energies, the baryon rapidity distribution exhibits a pile-up at 
mid-rapidity ||134| , |143|| (see figure |j). Whether the shape of the dN/dy distribution at 
SPS energies is flat or shows two bumps is currently not fully resolved (the SPS-data are 
preliminary). There are indications, however, that with rising beam energy the scaled 
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Figure 5. Comparison of baryon stopping in HIJING (solid), HIJING/B (dashed) 



and HIJING/B with "ropes" with various data [ 131 1 . 



dN / d(y I y p ) distribution stretches over the increasing rapidity gap between projectile 
and target; this can be seen in figure 0. Recently the NA49 collaboration ||144|| reported 
a A rapidity distribution which may be peaked strongly at mid-rapidity for Pb+Pb 
at 160 GeV/nucleon. This finding, however, is preliminary and in disagreement with 
equally preliminary results by the WA97 collaboration [|145|| , which indicate rapidity 
densities for A's lower by a factor of 2 — 3. 

Transverse energy measurements at the AGS [|146|| indicate that the transverse 
energy E T increases by 50% faster than predicted by an independent nucleon-nucleon 
interaction model when going from a light system (Si+Al) to a heavy system (Au+Au). 
In terms of a microscopic hadronic model this can be understood as a strong increase 
in baryonic density in the initial reaction phase and a corresponding large increase in 
the volume of high density matter |)9, |147|1 . 

The A(1232) abundance has been measured via tt + — p correlations at the AGS by 



the E814 and E877 collaborations ||148| , |149|| . The pion spectra can be decomposed into 
a thermal contribution and a contribution due to A— resonance decays. The A(1232)- 
to-nucleon ratio at freeze-out was determined to be ~ 35% for central silicon nucleus 
collisions. Hence, one can conclude that a large fraction of the system resides in hadron 
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Au+Au at SIS Au+Au at AGS 
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Figure 6. Time evolution of particle multiplicities (scaled with the number of 
incident nucleons) for central Au+Au collisions at 1 GeV/nucleon (SIS) and at 10.6 
GeV/nucleon (AGS). At SIS energies, only about 10% of the nucleons are excited to 
resonances whereas at AGS energies the degree of excitation exceeds 50%. For a time- 
span of up to 10 fm/c the baryons are in a state of A— matter. The figure has been 
taken from [137]. 
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Figure 7. Baryon stopping in central Pb+Pb collisions at 160 GeV/nucleon. 
Left: data by the NA49 collaboration (preliminary, figure taken from |l38| ) The solid 
line shows the rapidity distribution for net-baryons which can be decomposed into 
contributions from net-protons and net-A. For comparison the net-baryon distribution 
for central S+S collisions is also plotted (triangles). Right: UrQMD prediction 
compared to the same data (figure taken from |ll(| ). 



resonances, which produce most of the observed hadrons by their decay ("feeding"), 
after the resonances have decoupled. The dense state before this decay can therefore 
be called "resonance-matter". It exists due to the inertial confinement of energy and 
baryon number in the early phase of the reaction (see figure ||). 
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At CERN/SPS measurements of E T have been used to estimate the created energy- 
density: For 200 GeV/u S+Au central collisions ||150|| e reaches ~ 3 GeV/fm 3 . For the 
Pb+Pb experiment at 160 GeV/u similar values were extracted ||151| , but over a much 
larger volume. The reader is reminded here of the sensitivity of these extracted values 
on the hadron production time To, which is uncertain to at least a factor of two. 

Pb + Pb, NA49 Preliminary 
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Figure 8. Rapidity density of negative hadrons for central collisions. The circles 



represent the preliminary Pb+Pb measurement by the NA49 collaboration [153] 



whereas the squares are from the NA35 S+S experiment 152]. The latter are scaled 
by a factor of 6.6 , which corresponds to the relative number of participants. This 
scaled sulfur distribution agrees well with the lead distribution. The figure has been 



taken from [157] 



The mass- dependence of the rapidity distribution of produced particles, i.e. pions 
or kaons, can also used to search for scaling violations which could signal the onset of 
new physics phenomena. A comparison of the negative hadron rapidity distributions 



for S+S ||152|| with those for the lead on lead run ||153|| shows that the preliminary lead 
data can be matched by scaling the sulfur data with a factor 6.6 , close to the relative 
number of participant nucleons in central lead-lead collisions (A Pb /A s = 208/32 « 6.5) 
il53| (Figure |). 



Discussion 

The form of the measured baryon rapidity distributions shows experimentally that the 
central rapidity region up to Ei a b ~ 200 GeV/nucleon is not net-baryon free, in contrast 
to what had been expected in most early papers. Rather strong stopping as assumed 
first in hydrodynamic model studies P3"| , is observed. Therefore, results of theoretical 
analyses, which rely heavily on a net-baryon free mid-rapidity region with zero baryo- 



16 



chemical potential have to be taken with care. The quantitative measurements of the 
A-dependent stopping of baryons is one of the most important results of the AGS and 
SPS measurements. 

If the preliminary findings of a strongly peaked A rapidity distribution 



144 



and a rather broad Gaussian or flat baryon rapidity distribution ||153|| by the NA49 
collaboration are both confirmed, then this would be a hard obstacle for models which 
rely on global thermal equilibrium (plus flow) for the description of the final state of 
the the reaction pi, [[55], |5| . 



Simple "first collision models" without rescattering [|114j , |111|, |112| , |1 13| do not 
suffice to reproduce the data, whereas transport theory has correctly predicted the 
observed degree of baryon-stopping ||128| , |129| , |9~9|, |126| , |132| , |110|| . An alternative 
mechanism of baryon stopping based on diquark breaking [|158|1 , is also able to describe 
the corresponding experimental data, in contrast to the simple first collision approach. 
These models extrapolated to RHIC energies imply that even at y/s = 200 GeV/nucleon 
the dense matter is created with baryon density ~ 2p at r ~ 1 fm/c. In ||131|| the beam 
energy dependence of the initial baryon density is estimated to vary as 1/s 1 / 4 . 

The energy densities of e ~ 3 GeV/fm 3 estimated (with a factor of ~ 3 uncertainty) 
from rapidity distributions of produced particles indicate that part of the system may 
have entered the predicted state of deconfinement [E9 . Hadronic transport models, 



however, predict or reproduce the measured rapidity distributions, if baryon and meson 
rescattering and particle production via string decay ||128| , |129| , |116| , |9~9| , |110|| are included. 
Also hadronic models which include multi-quark droplets ||118|| above e cr it seem to give 
similar results. 

The inclusion of string excitations, collisions and decays are a first step towards 
modeling the parton/quark substructure of hadrons. In this sense these models go 
beyond what one would term purely hadronic model: Figure ^] shows y/s distributions 
of baryon-baryon interactions for Au+Au collisions at AGS and S+S collisions at SPS 
energies ||110|| . At the AGS, the collision spectrum is dominated by collisions of fully 
formed baryons. It exhibits a maximum at low energies, \/s ~ 3 GeV. Approximately 
20% of the collisions involve a diquark, i.e. a leading baryon originating from a string 
decay. In contrast to the heavy system at AGS, the collision spectrum for S+S at SPS 
exhibits two pronounced peaks. They are dominated by full BB collisions, one peak 
at the energy of initial projectile-target collisions, and one peak in the low ("thermal") 
energy range. Approximately 50% of the BB collisions, most of them represented by the 
bump at intermediate yfs values, involve diquark- or constituent quark-collisions with 
baryons. 

The linear scaling behavior of the mass dependence of the negative hadron rapidity 
distributions precludes a strong pQCD minijet component at these energies. We note 
that the agreement of the VNI parton cascade model with the Et systematics at SPS 
may be due to adjusting a strongly model dependent soft beam jet component to 



fit proton proton data. This issue is important because in ref. [ 160 1 it was claimed that 
at SPS already a partonic energy density > 5 GeV/fm 3 was created. The A 1 scaling of 
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BB collisions in Au(10.6GeV/c)Au 



BB collisions in S(200GeV/c)S 




Figure 9. -E^M distribution for baryon baryon collisions in central Au+Au reactions 
at the AGS (left) and in central S+S reactions at the SPS, calculated with the UrQMD 
transport model [11C ]. 



Et therefore constrains very strongly against hard scattering models used for example 
in [|161|| to argue for a QGP interpretation of J/ty suppression. We return to this point 
later. 

The experimental results demonstrate that highly excited dense matter is formed 
at mid-rapidity. They prove that a new state of elementary matter has been created. 
However, the inclusive central distributions do not give a clear and decisive answer to 
the question of whether this matter is predominantly of hadronic or quark nature. 

3.2. Creation of high temperatures: particle spectra 
Theoretical concepts 

The hot, dense reaction zone consists of slowed down incident nucleons and produced 
particles. The fireball model considers these hadrons as a mixture of ideal gases in 
thermodynamic equilibrium. For temperatures above 50 MeV and moderate densities, 
the Fermi- and Bose-Einstein-distribution functions for baryons and mesons (except for 
the pions) may be approximated by a Maxwell-Boltzmann distribution ||162| , |163|1 with 



the temperature T and the chemical potentials /Xj (connected to conserved quantum 
numbers i) as only free parameters. 

Kinetic equilibration is thought to be visible predominantly in the transverse 
degrees of freedom; therefore, transverse momentum or transverse mass distributions 
are used to extract temperatures from the spectral slopes. 

It has been suggested that abnormal nuclear matter, e.g. a QGP, may be observed 
via a secondary, high temperature component in the particle spectra or via a shoulder 



in the pion multiplicity distributions [164 



It has also been suggested that the equation of state, that is the energy density 
e vs. temperature T, can be probed experimentally by plotting the mean transverse 
momentum (p t ) vs. the rapidity density dN/dy or the transverse energy density 
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dN/dErp- If a phase transition occurs (i.e. a rapid change in the number of degrees 
of freedom) one expects a monotonously rising curve interrupted by a plateau: This 
plateau is is caused by the saturation of (p t ) during the mixed phase. After the phase 
transition from e.g. color singlet states to colored constituents has been completed ||165 



the mean transverse momentum rises again. However, detailed hydrodynamical studies 
166| , |167|| showed that the plateau is washed out due to collective flow. 

Collective (radial) flow [[54], |182| , |164|| as well as feeding from resonances strongly 



influence the shape of the particle spectra [|164|, [168| , |169| , |170| , |171| , |172| , |173| . For 
light composite particles, such as deuterons, the influence of collective flow is visible 
in a shoulder-arm shape of the transverse momentum spectra [p.64j| . This can be seen 



in figure [HJ. In order to account for flow effects, the spectra can be fitted with a 
thermal distribution including collective flow. The temperature T and the transverse 
flow velocity fit are the fit-parameters. The shapes of the velocity profile and density 
profile at freeze-out should enter as additional degrees of freedom in the analysis. Usually 
a box shaped density profile and a linearly increasing transverse velocity profile are 



assumed ||164| , |168| , |155| , |174|| . This results in severe distortions into the analysis, as 



discussed in the following [175 
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Figure 10. RQMD prediction of transverse mass spectra for protons and deuterons 
in central Au+Au collisions at the AGS compared with preliminary data by the E866 
collaboration. For deuterons a shoulder is visible in the low mt range of the spectrum. 



This structure is due to collective flow. The figure has been taken from [173] 



When extracting temperatures and flow velocities from microscopic calculations, 
the system is divided into cells and the local transverse and longitudinal velocity 
distributions are analyzed |125| , [L73| , |176| , |137|| . The temperatures extracted via a global 
two parameter fit are more than a factor of two higher than the temperatures gained from 
such a microscopic analysis at beam energies in the 100 MeV/nucleon to 10 GeV/nucleon 
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regime [ |176|| . The reason for this discrepancy lies mostly in the assumed shape of the 
freeze-out density profiles. 

Whereas a linearly increasing transverse freeze-out velocity profile seems a tolerable 
assumption, the shape of the freeze-out density profile has - due to collective flow - 
a Gaussian shape (centered at r t = 0), rather than the usually assumed box-shape 
distribution. When realistic density and velocity profiles are used, one finds that the 
high m t components of the particle spectra reflect contributions of large collective flow 
effects (i.e. the high expansion velocity). This analysis yields substantially lower values 
for the temperatures T. Such microscopic analyses of the spectra of protons, mesons 
and light composite particles at AGS energies show also that (3 t and T depend on the 
mass of the particle ||173| , |177| . 
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Figure 11. Excitation function of temperature T and average transverse expansion 



velocity f3 t . The figure has been taken from [178| 



Experimental status 

Data taken at the AGS with Si beams [140] seem on first sight to be consistent with 
an expanding, hadro-chemically and thermally equilibrated system with a temperature 
of 130 ± 10 MeV and a transverse flow velocity of (3 t « 0.36 |55|, [L74|. CERN SPS 
data with S beams have been fitted in the same fashion, with apparent temperatures 
around 150 MeV and flow velocities between 0.35 and 0.41 [170, |156|| . Figure [11] shows 
the extracted excitation function for the temperature T and the average transverse 
expansion velocity f3 t [ 178 |, including also SIS and BEVALAC data. 
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Figure 12. Mass and collision system dependence of the inverse slope parameter T 
measured by the NA44 collaboration [179] (left) and calculated by the UrQMD model 

(right) pa. 



In order to disentangle collective flow contributions from thermal motion, the 
dependence of the slope parameter T s i (which includes collective flow effects) on the 
collision system mass and the particle mass has been studied by the NA44 and NA49 
collaborations ||179| , |153| | at the SPS. Results can be seen in figure In proton-proton 
collisions obviously no collective effects are visible and an inverse slope parameter of 
T s i,pp = 145 MeV is extracted for all analyzed particle species (ir,K and p). When 
going to heavier collision systems, collective flow effects become obvious: the inverse 
slope parameter T s i increases with the mass of the emitted particle (see figure 12). 
Empirically one finds T s i = T s i tPP + m ■ {f3 t ) 2 . (3 t is the mean expansion velocity which 
depends on the mass of the collision system and m is the mass of the particle analyzed. 

The constant T s ^ pp in the empirical result therefore hints at the predicted limiting 



temperature | 181| , |164|| of 140 < T < 200 MeV. The observation of the increase of the 
flow effects for massive particles and heavy collision systems had been predicted with 
early hydrodynamical and microscopic calculations [|182j , |164| , f43[ , |183| , |187| , |180| . 



Discussion 

The data at AGS and CERN seem compatible with a hadro-chemically and thermally 
equilibrated system. However, this does not mean that the system necessarily evolved 
through thermal and chemical equilibrium states ||1 761, |184|, |185|. The fits to the spectra, 



with the temperature T and the transverse expansion velocity (5t as parameters, have 
to be performed with great care. There is a broad range of T and f3 t values which 
are compatible with the same spectrum [|164j , [170| , |171| , |176|| , where the temperature T 
depends crucially on the freeze-out density and velocity profiles - at least in the case of 
composite particles such as deuterons and tritons [|175| , [L73| , |186| , |180[ . 

The finding of one global freeze-out temperature T and velocity (3 t ||155| , |174| 



is to be contrasted with the independent analysis based on RQMD and UrQMD 
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calculations of spectra of light composite particles ||175| , |173| , |184j| and on spectra of 
mesons [ 187 , |184| , |110|| . These models are well able to reproduce the data and the analysis 
indicates different values (with a variation of ~20%) for (3 t and T, depending on the 
mass of the particle. The simplified thermal plus flow model ||164| , |170| , |155| , |174|| should 
not be taken literally. In reality we expect a complicated space-time dependent non- 
equilibrium freeze-out, details depending on inelastic production and absorption cross 
sections. In particular, the anti-baryon annihilation cross sections play an important 
role, as will be discussed below (section |3.3|) . Furthermore, flow of mesons vs. baryons 
188| . |189fl in opposite directions clearly indicate strong deviations from the single source 



fits as discussed in the following section. 

Recently, the WA98 collaboration reported ir° spectra in Pb+Pb reactions for p t up 
to 4 GeV/c |[L90|| . The data could be fit well by hydrodynamical models ||191| . However, 
it was found in ||192|| that the data were well reproduced by the QCD parton model. In 
this sense (the non-equilibrium) quark plasma is seen in the high p t spectra. However, 
as emphasized in |[L93||, at SPS the parton model is hypersensitive to models for soft 



multiple collisions. Hydrodynamics just happens to be one of the soft multiple collision 
models that can account for the data. 



3.3. Transverse collective radial and directed flow 
Theoretical concepts 

The excitation function of transverse collective flow is the earliest predicted signature 
for probing compressed nuclear matter [|57], |54 |. It has been shown that the excitation 
function of flow is sensitive to the EoS and can be used to search for abnormal matter 
states and phase transitions ||194| , |43| , |195| . 



In the fluid dynamical approach, the transverse collective flow is directly linked to 
the pressure of the matter in the reaction zone: 

With the pressure P(p, 5*) (depending on the density p and the entropy S), one can 
get a physical feeling for the generated collective transverse momentum p x by writing it 
as an integral of the pressure acting on a surface and over time [ |196|| : 



Px 



J J P(p, S) dAdt 



(2) 



where dA represents the surface element between the participant and spectator matters 
and the total pressure is the sum of the potential pressure and the kinetic pressure: The 
transverse collective flow depends directly on the equation of state, P(p, S). 

Directed collective flow has been predicted by nuclear fluid dynamics (NFD) 
54] , |197| , |198| , |199|| . Microscopic models such as VUU (Vlasov Uehling Uhlenbeck), 



and QMD (Quantum Molecular Dynamics) have predicted smaller flow than ideal 
NFD, these models show good agreement with viscous NFD [[L26|| and with the 
experimental findings [200, 201, 124, [202|| . It has been discovered initially at the 



the BEVALAC f20|, [204], |05[ for charged particles by the Plastic-Ball and Steamer 



Chamber collaborations p06| , at SATURNE by the DIOGENE collaboration |207| and 
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Figure 13. Excitation function of directed transverse flow. Left: prediction in the 
framework of nuclear hydrodynamics 219 , 216 1 , with and without dcconfinemcnt phase 
transition. In the case of a phase transition a minimum in the excitation function is 
clearly visible. Right: Data compilation taken from [443] 



has been studied extensively at GSI by the FOPI [£0|, [I7§, LAND P09| , TAPS |gTD 
and KaoS [2TT| collaborations. 



One has to distinguish two different signatures of directed collective flow: 



a) The bounce-off [ |f 97|1 of compressed matter in the reaction plane and 

b) the squeeze-out [|198|1 of the participant matter out of the reaction plane. 

The most strongly stopped, compressed matter around mid-rapidity is seen directly in 
the squeeze-out ||212|| . A strong dependence of these collective effects on the nuclear 
equation of state is predicted 11202}] . For higher beam energies, however, projectile and 
target spectator decouple quickly from the reaction zone, giving way to a preferential 
emission of matter in the reaction plane, even at mid-rapidity [[213|| . An excitation 
function of the squeeze-out at midrapidity, possibly showing the transition from out of 
plane enhancement to preferential in-plane emission has been predicted to enhance the 
sensitivity to the nuclear equation of state ||214| , |110|| . 

Apart from the above discussed directed flow, the so-called "radial" , i.e. undirected, 
flow component can be used for simplicity (azimuthal symmetry) (Tl 



161. It has to 



be taken into account for the interpretation of particle spectra used for temperature 
extraction which may drop by as much as a factor of 2. 

Due to it's direct dependence on the EoS, P(p, T), flow excitation functions can 
provide unique information about phase transitions: The formation of abnormal nuclear 
matter, e.g., yields a reduction of the collective flow ||194| . A directed flow excitation 
function as signature of the phase transition into the QGP has been proposed by several 



authors [|43], |119|1 . A microscopic analysis showed that the existence of a first order phase 



transition can show up as a reduction in the directed transverse flow ]p!2|| . 

For first order phase transitions, the pressure remains constant in the region of the 
phase coexistence. This results in a vanishing velocity of sound c s = J dp/de. 



23 



Pb(160AGeV) Pb, b=0-3fm 




2 4 

Rapidity 



10 



_^ 4 

2 

7 ° 

O -2 
-4 

-6 



-10 



Protons J 

r ■ 




V 


-A : 


Y 








) 1 2 




CO 


4 5 ( 



rapidity 



Figure 14. RQMD 1.05 prediction of collective sideward flow for the system Pb+Pb 



at 160 GeV (left, figure taken from [ 129 1 ) . The rhs shows data on directed and elliptic 
flow vs. rapidity by the NA49 collaboration (figure taken from |j23l[ ). 



The expansion of the system is driven by the pressure gradients, therefore expansion 
depends crucially on c 2 s . Matter in the mixed phase expands less rapidly than a hadron 
gas or a QGP at the same energy density and entropy. In case of rapid changes in the 
EoS without phase transition, the pressure gradients are finite, but still smaller than for 
an ideal gas EoS, and therefore the system expands more slowly [ |166| , |167| . 

This reduction of c 2 s in the transition region is commonly referred to as softening 
of the EoS. The respective region of energy densities has been called the soft region 
|183| , |215| , |216| , [217 . Here the flow will temporarily slow down (or possibly even stall). 
Consequently a time delay is expected in the expansion of the system. This prevents 
the deflection of spectator matter (the bounce-off) and, therefore, causes a reduction of 
the directed transverse flow |[218| , |219|| in semi-peripheral collisions. The softening of the 
EoS should be observable in the excitation function of the transverse directed flow of 
baryons (see figure |T3|). 

The overall decrease of p x seen in Fig. [13| for > 10 GeV both for the hadronic 
and the QGP equation of state demonstrates that faster spectators are less easily 
deflected (because A and t in equation ^ are decreasing with Ei a b) by the hot, expanding 
participant matter. For the QGP equation of state, however, these one-fluid calculations 
show a local minimum in the excitation function, at about 6 GeV/nucleon. This can be 
related to the QGP phase transition, i.e. to the existence of the soft region in the EoS. 

The limitation of one-fluid hydrodynamic calculations is that they assume 
instantaneous thermalization. This becomes unrealistic for increasing beam energies 
since due to the average rapidity loss of only one unit per proton-proton collision, 
nucleons require several collisions for thermalization. A more realistic three-fluid 
calculation, in which the third fluid represents the fireball of produced particles and 
only local thermal equilibrium is assumed, yields much lower flow values - even without 
a first order phase-transition ||120|| . The position of the minimum (the magnitude of 
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the overall effect) therefore strongly depends on the degree of stopping (i.e. which type 
of fluid-dynamical model is employed) and on the details of the chosen EoS and phase 
transition parameters. 

Taking the finite volume of the reaction zone into account, one finds that 
fluctuations hinder a sharp separation between the QGP-phase and the hadronic phase 
and lead to a rounding of the phase transition ||84|| . For realistic reaction volumes the 
softening of the equation of state is reduced considerably and thus the minimum-signal 
in the flow excitation function is washed out. 

A second order phase transition may not exhibit this minimum in the flow excitation 
function: The existence of a minimum in p Xy dir{Ei a b) is rather a qualitative signal for 
a strong first order transition. If such a drop of p x ,dir{Ei a b) is observed, it remains 
to be seen which phase transition caused this behavior: a hadron-quark-gluon phase 
transition or, e.g., a resonance matter - ground state matter phase transition in confined 
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Figure 15. Excitation function of directed transverse flow for different collision 



systems from fOO MeV/nucleon up to 2 GeV/nucleon. Left: VUU prediction [101 
for the system Au+Au. Right: data from the Streamer Chamber, Plastic Ball 
and EOS experiments at the BEVALAC and from the FOPI experiment at SIS 
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178| (the figure has been taken from [178|). The data has been 



A| ) to account for the different collision systems. 



Experimental status 

Collective flow measurements have first been performed at the BEVALAC [ 203 , |204| , 205 



for charged particles by the Plastic-Ball and Streamer Chamber collaborations. A more 
detailed investigation of the excitation function between 0.1 to 1.2 GeV/nucleon for 
Au+Au has been performed by the FOPI, KaoS, LAND and TAPS collaborations at 
GSI |7|, |209], gig, and the EOS-TPC collaboration at LBNL |22| (see figure [15]). 
At 10.6 GeV collective flow has recently been discovered by the E877 collaboration 

d((E x ) I (Et)) /dr] for different centrality bins. 



222, 223 



Figure [16] shows dvi/dr] 
The E895 group has measured the flow excitation function for Au+Au at the AGS in 
the energy range between 2.0 and 10.6 GeV/nucleon ||224|| . Their data show a smooth 
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Figure 16. Transverse collective flow measured at the AGS by the E877 
collaboration ]223[ . Plotted is dvi/drj — d((E x ) / (Et)) / drj which is a similar quantity 
as d((j> x ) j (pt))dy, used at lower beam energies. 



decrease in (p x ) from 2 to 8 GeV/nucleon and are corroborated by measurements of the 
E917 collaboration at 8 and 10.6 GeV/nucleon |225| . 



The E895 collaboration has also measured an elliptic flow excitation function 
indicating a transition from out-of-plane enhancement (i.e. squeeze-out) to in-plane 
enhancement around 5 GeV/nucleon ||226 . 

At CERN/SPS, the existence of undirected flow has been deduced from a combined 
analysis of particle spectra | 227 , 156 | and HBT correlations ||228|| (see also sections ET2" 
and 



First observations of a directed transverse flow component have been reported by 
the WA98 collaboration | 22H , 230 using the Plastic Ball detector located at target 
rapidity for event plane reconstruction. They show a strong directed flow signal for 
protons and "antifiow" for pions, both enhanced for particles with high transverse 
momenta. The same findings have been reported from the NA49 collaboration, which 
due to its larger acceptance allows for an even more detailed investigation. They 
report a quite strong elliptic flow signal near mid-rapidity at 160 GeV/nucleon |[231| 



(see figure 14). 



Discussion 

An observation of the predicted local minimum in the excitation function of the directed 
transverse flow [[21 9| , |216|| would be an important discovery, and an unambiguous signal 
for a strong phase transition in dense matter. It's experimental measurement would 
serve as strong evidence for a QGP and a strongly first order deconfinement transition 
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at non-zero baryon density. 

A strong experimental effort at the AGS and SPS has led to the discovery of 
flow even at these ultra-relativistic energies. The search for the minimum-signal in 
the excitation function is under way. 

The absolute values for the ideal NFD prediction of directed flow shown in figure 
overestimate the experimental values considerably |[21G| , |217|| due to lack of viscosity 
|126| . The position of the minimum in p X; dir(Ei ab ) depends on the EoS - therefore it 



is by no means clear where (in Ei a b) the deconfinement phase transition will occur. 
Furthermore, finite volume corrections reduce the softening of the equation of state and 
might reduce the minimum-signal considerably H84| . 

The combined efforts of the FOPI and EOS/E895 collaborations will allow to map 
experimentally the region from 0.1 GeV/nucleon to 10 GeV/nucleon. However, the 
current data show a smooth decrease in the flow from 2 to 10 GeV. This seems to 
favor a hadronic scenario without a phase transition. An experimental search for this 
outset of flow in the energy range between 10 GeV/nucleon and 200 GeV/nucleon seems 
necessary. 

The recent measurement of the squeeze-out excitation function between 2 and 8 
GeV/nucleon may offer a new approach for studying the nuclear equation of state 

The comparison of proton spectra with 0-meson spectra may help to disentangle 
"early", QGP-related flow components from "late", hadronic contributions. Transport 
model calculations have shown that the 0-meson decouples much earlier from the system 
12 fm/c) than the nucleons [232|. Since both particles have approximately the same 



mass, their "thermal" motion and undirected flow components should be identical and 
any differences in the spectra should arise only through the additional interaction the 
nucleons suffer in the later reaction stages [[23 2j . 



3.4- Space time pictures of the reaction: HBT source radii 
Theoretical concepts 

Intensity interferometry of identical particle pairs, such as tttt, KK or pp pairs, can 
be used to extract information about the space-time dynamics, freeze-out volume 
and reaction geometry of heavy-ion collisions. The method was originally devised 
by Hanbury-Brown and Twiss to measure the angular diameter of a star using the 
correlation of two photons ||233|1 : 



The probability of detecting two photons in coincidence in two different detectors 
is correlated to the relative separation between the two detectors. This correlation is 
connected to the angular diameter of the emitting source. This effect is commonly 
known as the Hanbury-Brown- Twiss (HBT) effect. It has been also observed in proton- 



antiproton annihilations 234 



By applying the HBT-measurements to particles emitted in heavy ion reactions, 
such as protons, pions or kaons, the two particle correlation function yields the 
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longitudinal and transverse radii as well as the lifetime and flow pattern of the emitting 
source at the moment of freeze out [ |235| , |236| , |237| , |238| , [239| , [240|| . The inverse widths 



R out of the "out" correlation function and R S ide of the "side" correlation function can 
be used to extract a measure for the duration of particle emission (R 2 out — R 2 side ) and the 
transverse size of the source (R S ide) [ |238j , |239|| . 



The prolonged life-time of the collision system in the mixed phase, which has 
already been discussed in section |3.3j , can be observed through an enhancement of 
the ratio of inverse widths (R ou t/ Rside) of the two particle correlation function in out- 
and side-direction [|38|, [Ml], ^42), |17| For energy densities estimated to be reached 



in Pb+Pb collisions at the CERN/SPS one expects R out /R S ide ~ 1.5 — 2 [[21711 



Inclusion of the decays of long-lived resonances may however reduce the R ou t/ Rside 

ratio m, m, m, m nzi- 



Final state interactions between non-identical particles can provide information not 
only about the duration of the emission but also about its time-ordering. It has recently 
been shown that an anisotropy in the space-time distribution of emitted particles reflects 
in the directional dependence of unlike-particle correlations (e.g. p — K) and can 
thus directly be used to measure the sequence of the emission of particles of different 



types ||248|| . Applying this technique to the correlation between a strange and an anti- 
strange particle (e.g. K + K~ interferometry) [|249|| may result in the direct observation 
of the strangeness distillation process ||250|| (see section pT6| ) . That process - which 
is instrumental to the formation of so-called strangelets - predicts an enrichment of 
s quarks in the quark phase while the s quarks drift into the hadronic phase. The 
resulting time-ordering of the freeze-out for strange and anti-strange particles is to be 
compared to the (different) emission times due to the different mean free paths in a 
purely hadronic scenario | |249[ | . 

A combined analysis of single- and two-particle spectra can yield a rather complex 
reconstruction of the geometry and dynamical state of the source at freeze-out [[251| . 
This information can be used as a powerful test for dynamical simulations of the collision 
process. 

Experimental status 

K + K + and K~~K~ measurements at the SPS ||253|| show similar radii around 2.7 (±10%) 



fm for the system S+Pb. Since the i^~-nucleon interaction cross section is far larger 
than the i^ + -nucleon cross section, this result indicates that the dominant interaction 
for kaons in the later reaction stages (close to freeze-out) at SPS energies are K-ir 
interactions ||255|| . At AGS energies, the situation might be different: the baryon to 
meson multiplicity ratio is approximately one, there. A detailed analysis has yet to be 
performed. Radii extracted from nir correlations are larger than those from KK, both 
at AGS and SPS energies [|253| , |256|1 . The differences are caused by different interaction 



cross sections and resonance decays [|257| , [258| , plus the effect of collective expansion 



240, 259 . More theoretical work is needed to separate these effects. 



For central collisions of heavy systems the extracted transverse radii are on the 
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Figure 17. Systematics of HBT radius parameters, compiled from data by the NA44. 



The figure has been taken from [261 



order of 5 to 7 fm for the vrvr (both, AGS and SPS) and 3 fm (AGS) to 4 fm (SPS) for 
the KK system H, [261], [262], [263], |6§ . 



The longitudinal and transverse radii measured at the SPS are larger than the 
respective radii of both, the projectile and the target, indicating an expansion of the 
system prior to freeze-out [[26 5| , |251|| . Data with the sulfur beam at the SPS show that 
the longitudinal radii measured as a function of rapidity |[266|| could be fitted by a boost 
invariant longitudinal expansion [267|. Recent data taken with the lead beam by the 
NA49 collaboration confirm this finding for the Pb+Pb system ||262| , [268| , |264| , |228|| . 

Transverse source radii for tttt show a decrease from 4 to 2.5 fm for S+S and from 
6 to 4 fm for Pb+Pb, respectively, with increasing transverse momentum of the pions 
(see also Figure [17|) [ (266| , |262| , [263| , |264|| ; this behavior is to be expected in the presence 
of transverse flow |237| , |238| , |239| . Alternatively, it can also be explained by microscopic 
models which predict that high p t particles are emitted in the early reaction stages (by 
heavy resonances or strings) and low p t particles (which have rescattered more often) 
have late freeze-out times ||269| , |137|| . Measurements of R out /R S ide indicate values on the 
order of 1 [ |265|| . For the Pb+Pb system, a duration of emission of about 3-4 fm/c has 
been reported with the lifetime of the source being r « 8 fm/c ||26 1 . 



Source radii cannot only be studied as a function of transverse momentum or beam 
energy, but also as a function of impact parameter related quantities, such as the number 
of participant nucleons. The latter analysis can be performed either by comparing 
central events of different systems or by comparing different centralities in very heavy 
systems. 

Exciting preliminary results have been reported at QM '96 for the system Au+Au 
at 10.6 GeV/nucleon by the E866 collaboration |[270|| , showing a dramatic increase of 
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40% in the source radius over the last 7% of highest centrality. More recent data-sets 
of the same collaboration suggest a more gradual increase with centrality |[271|| . The 
limited statistics, however, do not permit a final assessment, yet. 



Discussion 

At AGS energies hadronic transport models are well able to reproduce the measured 
source radii [|256j |260| - at SPS energies a full analysis has not yet been performed, but 



early comparisons showed at least qualitative agreement ||256| , |253 
A strong first order phase transition 



241|| and even an infinite order but rapid 



cross over transition [[21711 should result in a lower pressure, slower expansion and perhaps 
a long-lived evaporating droplet of QGP. The rather short lifetime of r 8 fm/c ||264|| for 
Pb+Pb at CERN/SPS suggests either the non-existence of such a low-pressure system 
or perhaps that the initial energy-density that is needed to create a QGP is much higher 

Sim- 

HBT-interferometry shows thusfar no evidence for the characteristic time delay of 
QGP formation up to SPS energies. The main complication of HBT analysis in nuclear 
collisions is the existence of strong collective flow that (Doppler) distorts the interference 
pattern. The p t dependence of the HBT radii has become a useful tool to probe this 
aspect of the reaction dynamics. It will be important to search for time signatures at 
RHIC and LHC. 



3.5. Remnants of hadronization: strangeness enhancement 
Theoretical concepts 

In proton proton collisions, the production of particles containing strange quarks is 
strongly suppressed as compared to the production of particles with u and d quarks 
|272| , [273j| . It has been argued that this suppression is due to the higher mass of the 
ss quark pair. The suppression increases with the strangeness content of the particles 
produced in proton proton collisions. 

In the case of QGP formation, ss pairs can either be produced via the interactions 
of two gluons or of qq pairs. Leading order a s pQCD calculations suggest that the 
second process dominates only for y/s < 0.6 GeV |[274| |. The time-scale of chemical 
equilibration of (anti-) strangeness due to gluon gluon interaction is estimated - also 
based on first order pQCD calculations - to be about 3 to 6 fm/c, depending on the 
temperature of the plasma ||275 |. 

Following this line of argument, the yield of strange and multi-strange mesons and 
(anti-) baryons has been predicted to be strongly enhanced in the presence of a QGP as 



compared to a purely hadronic scenario at the same temperature ||276| , |277 |. However, 
the estimated equilibration times may not be sufficiently rapid to cause a saturation in 
the production of strange hadrons before QGP freeze-out. 

In particular, assuming low chemical potentials, fi^ ~ fi u ~ = fi s and a 
temperature T higher than the strange quark mass m s , the densities of all quarks 
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and anti-quarks are nearly the same in the QGP. Hence, the probability of forming 
anti-hyperons by combining u, d and s quarks is nearly the same as the probability of 
forming strange and non- strange baryons by combining u, d and s quarks if the freeze 
out process is rapid and annihilation can be neglected. 

In contrast, the production of an antihyperon-hyperon pair produced in nucleon 
nucleon collisions is greatly suppressed by the Schwinger factor P78L |279 



since it is 

necessary to tunnel the massive diquark and the strange quark through the potential 
wall in the chromo-electric field with the string tension k « 1 GeV/fm ||280|| . The 
enhanced production of anti-hyperons (A, S, S and Cl) can therefore be used as a QGP 
signal in the case of zero chemical potential ||275|| . 

If a QGP is created in heavy ion collisions at AGS or SPS energies, it will most 
likely be characterized by nonzero chemical potentials /i u and fid- This results in the 
densities of u and d quarks being larger than those of the s and s quarks, which in turn 
are larger than the u and d densities. Due to these different abundances the s quark is 
more likely to combine with aword quark to form a K + or K° (or with two non-strange 
quarks to form a A or S, respectively) than it is for the s quark to recombine with a u 
or d quark thus forming a K°s or K~s. Therefore, in the QGP case the K + /tt + ratio 



in a relativistic heavy ion collision is different from the K jix ratio [[281 



The relative abundances of various strange particle species have been used for the 
determination of relative strangeness equilibration. To account for incomplete chemical 
equilibration, a strangeness fugacity j s is introduced in a thermo-chemical approach 
||282| , [49| , |283| , |284|| . One has also compared the measured ratios and the connected 
thermodynamic variables (such as T, \i B and the entropy) with calculations, either 
assuming a hadron gas scenario or a QGP scenario including some hadronization scheme 

There are certain drawbacks to the line of argument presented above: The strange 
particle abundances, after freeze out from a QGP, are very close to those of a fully 
equilibrated hadron gas at the same entropy content ||287|| . The reason is ||281| , [288] , [289| 



that the volume of a hadron gas of the same total energy has to be larger due to the 
smaller number of available degrees of freedom. Consequently, one must expect that the 
abundance of strange quarks is diluted during the hadronization process. This dilution 

where gluons hadronize by 



effect is clearly seen in hadronization models ||277| , p29C 
conversion into quark-anti-quark pairs, which predominantly feed the final pion channel. 
As a consequence, the K/ir ratios are significantly reduced. 

Furthermore, the computation of particle abundances in the QGP and the hadron 
gas scenario are mostly based on the assumption of chemical and thermal equilibrium 
(a non-equilibrium calculation has been published in |[290||). For the hadronic case these 



assumptions cannot be justified: It has been shown via rate equations [ |277| , |281|| that 
the strangeness equilibration time exceeds the reaction time of a heavy ion collision by 
at least one order of magnitude. 

Strangeness production in the hadronic scenario is a non-equilibrium process. In 
the early (pre-equilibrium) reaction stages, typical longitudinal momenta are much 



31 



higher than in the case of a thermal momentum distribution. This leads to enhanced 
strangeness production ||291| , |232|| . The system then cools down in the course of the 
reaction. It's final "equilibrium" temperature is therefore only partly connected to the 
measured strange particle yields and spectra. 



Experimental status 

An enhancement of the K/tt ratio has been measured both at the AGS and at the SPS 
292| . At the AGS, K + /ir+ w 0.2 and K-/ir~ w 0.04. Furthermore K + /K~, A/A 

At the SPS, 
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and p/p production ratios have been measured at the AGS 
enhanced production of (anti)hyperons, such as A, S, Q and Q, has been observed and 
ratios of A/A,S/S,S/A and S/A have been analyzed by the NA36,WA85 and WA97 
collaborations |9|, |95|, |9|, |29^, |9|, |9|, |300| . The WA94 collaboration has 
measured antihyperon ratios (i.e. the S/A ratio) in pp, pA and AA reactions. They 
find a smooth increase in the S/A value from pp over pA to AA ||301|| . 

Recently, very interesting values have been quoted for the A/p ratio. It has been 



measured by the NA35, NA49, E866, E878 and E864 collaborations |03], [153|, |04], [305 



Since it only contains newly produced anti-quarks, it may therefore represent a rather 
clean measure for the s/u quark ratio in the hot and dense matter. For pp and pA 
collisions this ratio is below 0.4, whereas in AA collisions preliminary analysis give 
values between 3 and 5 | |305| , |144| | - these values are so high that they could not be 
obtained in either a hadron gas or QGP model with reasonable values for T, fis and fis- 
The observed strong enhancement of multistrange (anti)hyperons (S, Q, S, and 
Q) from light to heavy collision systems at the CERN/SPS |Ul| |00], |02| surely 
constitutes on the experimental side the most intriguing evidence for a possible non- 
hadronic enhancement of strangeness. 



Discussion 

Hadronic models for particle production ||306| , |287| , |291| , |307| , |99| work quite well in 
the case of the observed K + /tt + enhancement |[292j 1 at the AGS (silicon beam). The 
reason for strangeness enhancement in a hadronic scenario is multistep excitation of 
heavy baryon and meson resonance states ||291|| . The AGS value of K + /n + m 0.2 is 



compatible with a strangeness equilibrated hadron gas ||174 . 

AGS data of K + /K~,A/A and p/p ratios can be fitted with an equilibrated 
hadronic fireball with fi s /T = 0.54±0.11 and fx B /T = 3.9±0.3 mil Gil . However, 



this does not mean that the system has always been in the hadronic phase, since an 
equilibrium state has no memory on how it has been produced. The system might 
as well have originated in the quark phase and then evolve along the phase-boundary, 
thereby hadronizing with varying combinations of T, /i# and The point is that these 
ratios provide actually very little information about the properties of the early time 
dense system. 

The ratios A/A, S/S, S/A and S/A measured by the WA85 and WA97 
collaborations ||295| , |296| , |297|| at the SPS can be fit in analogy to the ratios at the 
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AGS by an equilibrium hadron gas model with 7 S = 0.7, /i# = 0.24 and T = 180 MeV 
154j| . Besides the three parameters T, fig and us which are used in the grand-canonical 



formalism of statistical mechanics, the additional parameter 7 S accounts for incomplete 
saturation of strange particles in phase space. However, data can also be fitted with 



a hadron gas model and 7 S ~ 1 with fJ, s /T = 0.24 — 0.28 and [ib/T = 1.05 |[L56|| , 
respectively. The very same data can also be fit by an instantaneously hadronizing 
non-equilibrated QGP with strangeness neutrality and strangeness saturation 7 S > 0.7 

mum. 
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Figure 18. Left: final particle ratios computed in UrQMD (full circles) [316] and 
a non-equilibrium hadronization scenario (crosses) with initial conditions Ag 1 * = 100, 
Sy^init = 45; yinit = and bag constant £i/4 = 235 MeV |ll. The data (open 



circles) are taken from various experiments as compiled in [156]. Right: corresponding 
particle production rates as a function of time. Strong differences in the time-evolution 
of various particle ratios are observed. 



Are the extracted temperatures and chemical potentials really reliable in view of 
the simple, static, thermal ansatz? Even hadron production in high energy pp and 
pp collisions has been calculated by assuming thermal and chemical equilibrium and 
fits the data well ||311|| . The fit temperature lies around 130 MeV - 170 MeV, nearly 
independent of the center of mass energy of the incident particles. A 7 S value of ~ 0.5 is 
needed for the fit, indicating incomplete strangeness saturation already at the pp level. 
Does such a model make sense? The success of the fit can be interpreted as hadron 
production in elementary high energy collisions being dominated by phase space rather 
than by microscopic dynamics. 

The extrapolation of this conclusion to heavy-ion collisions, however, may not 
be valid: Even simple dynamical hadronization schemes [|312j| , where thermodynamic 
equilibrium between a quark blob and the hadron layer is imposed, reveal a more 
complex picture (see figure 0). Particle ratios can be reproduced nicely with the 
same number of parameters as in the static ansatz of a hadron gas in equilibrium, 
while the space-time evolution of the system shows strong changes of the strange and 
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baryo-chemical potentials due to baryon- and strangeness-distillery [|313| , |314|1 . Taking 
(boost-invariant) longitudinal hydrodynamical expansion into account, the interplay of 
the evaporation process and the hydrodynamical expansion (and vice versa) leads to 
considerably shorter lifetimes of the mixed phase as compared to scenarios without 
hydrodynamical expansion ||315| . It is very questionable, whether final particle yields 
reflect the actual thermodynamic properties of the system at any one stage of the 
evolution. 

Microscopic transport model calculations are in good agreement with the measured 
hadron ratios of the system S+Au at CERN/SPS [|316| , |317|1 . They show, however, that 
those ratios exhibit a strong rapidity dependence. Thus, thermal model fits to data may 
be distorted due to varying experimental acceptances for individual ratios. A thermal 
model fit to S+Au ratios calculated with the microscopic UrQMD transport model 
(and extracted within the same range of rapidity for all ratios) yields a temperature of 
T=145 MeV and a chemical potential of \lr = 165 MeV ||316|| . Hadron ratios for the 
system Pb+Pb are predicted and can be fitted by a thermal model with T=140 MeV 
and fiB — 210 MeV. Similar results have been obtained with RQMD |ftl7| . Analyzing 
the results of a non-equilibrium transport model calculations in the framework of an 
equilibrium model may, however, not seem meaningful. 

The smooth increase of the H/A ratio from pp via pA to AA reactions suggests 
that production volume and the degree of thermalization may not be relevant for the 
production of antihyperons. Already two overlapping strings (as typically produced 
in p+S reactions) are sufficient to yield strong deviations from the behavior in pp 
232|, |318|]. Clearly the detailed study of pA reactions yields important information 



on the production processes of antihyperons. 

However, so far all models fail to describe the recently reported, unusually high 
A/p ratio of ~ 3 — 5 (proton-proton collisions yield a ratio of 0.2 - 0.3). One 
possible explanation could be that the A have a far lower annihilation cross section 
than the p. This difference in the annihilation cross section might account for the 
dramatic A enhancement. A straightforward way to test this hypothesis would be the 
measurement of p and A (anti-)flow. For p's, a strong anti-correlation with regard to the 
"conventional" baryon flow is predicted ||189|| . This is due to their large annihilation cross 
section in dense matter. The same would only hold true for the A, if its annihilation 
cross section is correspondingly large. If, however, A anti-flow is not observed, this 
would serve as clear indication for a low A annihilation cross section. Thus the above 
explanation for the A enhancement as being due to smaller a ann ||319|| would be supported 
by independent evidence. 

Alternatively, the A/p enhancement could be explained by different medium 
modifications to the masses of non-strange and strange baryons which affect the 
production probabilities. If an attractive strange scalar condensate lowers the mass 
of the A in hot and dense hadronic matter, even below that of the p, this could account 
for the A/p enhancement [ |319] , 320 . 



Thermal model analyses assume constant freeze-out temperatures and chemical 
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potentials, but at least the more careful ones do not assume a static source; instead they 
allow for collective expansion flow. While the flow does not matter for an analysis of 
47r yields, it becomes indeed important when comparing the model to data from limited 
windows in momentum space. Unfortunately, no conclusion is possible unless the freeze- 
out surface is known. Most people use too simplistic isochronous (t =const.) freeze-out 
prescriptions which in fact correspond to a volume freeze-out. However, it has been 
shown in the framework of an expanding hadron-gas that freeze-out is not a state but 
a reaction stage and that the various equilibria (i.e. chemical and thermal equilibrium) 
necessarily break down in the final stages ||321|| . Microscopic model calculations support 



this picture of a complicated sequential freeze-out depending on reaction rates and 
particle species ||187| , |185| , |232|| : Even if some particles are in thermal and chemical 



equilibrium during the final stages of the reaction, the problem of how to disentangle 
the thermal contribution from the early pre-equilibrium emission would remain. This 
problem has not been addressed satisfactorily so far (see section |3~2|) . 

Hadronic transport models, which are based on a non-equilibrium scenario, however, 
are only able to describe the CERN/SPS (anti-)hyperon data by invoking non-hadronic 
scenarios such as color-ropes ||116|| , breaking of multiple-strings ||117|| or decay of multi- 



quark droplets [|118|| . Therefore, the (anti-) (strange-) baryon sector remains a topic 
of great interest. Specifically the strong enhancement of multistrange (anti)hyperons 
(S, Q, S, and Q) heavy collision systems such as Pb+Pb at the CERN/SPS is of great 
importance [ |300|| since it offers currently the best opportunity to discriminate hadronic 
from deconfinement scenarios in the sector of strangeness enhancement. 

3.6. Ashes of the plasma: strangelets and hypermatter 
Theoretical concepts 

The observed abundant production of strange baryons at AGS and SPS energies led 
people to speculate about implications for hypermatter (multi-hyperon clusters or 
strange quark droplets ) formation [ |322j , |323| , |324| , |325| , [326|| . Speculations about the 



existence of such multi-strange objects, with baryon numbers B > 100, have been 
around for decades, in particular within astrophysics. Such states are allowed for by the 
standard model, although so far their existence has not been proven in nature, e.g. in 
the form of strange neutron stars. 

Quark matter systems with A > 1 are unstable, if they only consist of u and d 
quarks, due to the large Fermi energy of these non-strange quarks. The system's energy 
may be lowered by converting some of the u and d quarks into s quarks (i.e. introducing 
a new degree of freedom). The energy gain may over- compensate the high mass of the 
s quarks - thus such strange quark matter (SQM) may be absolutely stable ||325| . 



Hadrons with B > 1 and S < have been considered even before the advent of QCD 
322 , 323 1 . However, first the development of the MIT Bag Model [ [327 ] allowed to model 



such states. Long hypermatter lifetimes (for hundreds of quarks and a strangeness per 
baryon ratio in the order of 1) have been predicted, up to 10 -4 seconds [ 324 |. Further 
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detailed investigations of small pieces of strange quark matter, so called strangelets, 
reveal possible (meta) stability for B > 6 [ |325| , |326| . The simplest strangelet is the 



H— dibaryon with zero charge, B = 2 and S = —2, which consists of 2u, 2d and 2s 
quarks, followed by the strange quark-a with 6u, 6d and 6s quarks [ |325| , |328| . 

For a QGP - hadron fluid first order phase transition with nonzero baryo-chemical 
potential, a mechanism analogous to associated kaon-pro duct ion yields an enriched 
population of s quarks in the quark-gluon phase, while the s quarks drift into the hadron 
phase [|250| , |329| . This strangeness separation results in the distillation of metastable 



strangelets only if the Bag constants are very small, B < 180 MeV/fm 3 [250] . 



Experimentally strangelets are distinguishable from normal nuclei due to their very 
small or even negative charge to mass ratio. The most interesting candidates for long- 
lived strangelets are lying in a valley of stability which starts at the quark-a and 
continues by adding one unit of negative charge, i.e. (A,Z) = (8,-2),(9,-3). . . ||330|| . Recent 
calculations indicate that positively charged strangelets seem only to exist for A > 12 
and very low bag parameters |[330|| . 

There exist, however, other forms of hypermatter with similar properties as 
strangelets: hyperclusters or MEMO's (metastable exotic multihyperon objects) consist 
of multiple A, X and S hyperons ||331| , and - possibly - also nucleons. The double- 



A hypernucleus \ A He has been observed long ago [ |332|| . Properties of MEMOs have 
been estimated using the Relativistic Mean Field model. MEMOs can contain multiple 
negatively charged hyperons, therefore they may also have zero or negative charge-to- 
mass-ratios. 

MEMOs or hyperclusters could form a doorway state to strangelet production, or 
vice versa: MEMOs may coalescence in the high multiplicity region of the reaction. If 
strangelets are stronger bound than "conventional" confined MEMOs, the latter may 
transform into strangelets. The cross sections for production of MEMOs in relativistic 
heavy ion collisions rely heavily on model parameters (e.g. in the in the coalescence 
model po and r ). The predicted yields are typically < 10 -8 per event [|331| , |333|| . 



Experimental status 

Strangelet searches are underway at the AGS ||334j , |335| , |336|| and SPS ||337| , |338| , |339| , [340 



|341| . So far no long lived (r > 10~ 7 s) strangelets have been unambiguously identified 
- the upper limits for the production cross sections established by the experiment are 
still consistent with theoretical predictions for short lived MEMOs since they cannot 
be tested in the present long flight path experiments. There has been a report of one 



candidate with Z = -1, N/Z = 7.4 GeV and r > 85 (is |41|, |3|, |5I2] . Therefore this 
exciting topic awaits more experimental effort. 

Current experiments are designed to detect strangelets with a small charge-to-mass 
ratio and rather long lifetime (r > 12 fis in the case of ||339| , |340| ) . The present 



experimental setups are hardly sensitive to the most promising long-lived and negatively 
charged strangelet candidates beyond the strange quark-a. Unfortunately, plans for 
extending experiment E864 at the AGS to look for highly charged strangelets with 
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B > 10 ||343|| cannot be followed because the AGS fixed target heavy ion program has 
been put to rest. 

Future experiments at collider energies (STAR at RHIC and ALICE at LHC) will 
be sensitive for short-lived metastable hypermatter, too ||504|, [344], |345| . 



Discussion 

Due to the possibility of creating MEMOs in a hadronic scenario and their possible 
subsequent transformation into strangelets, the formation of a QGP is not a necessary 
prerequisite for the creation of strangelets. The discovery of strangelets would therefore 
be no hard proof for a deconfmement phase transition. So far there seems to be no 
clear way to distinguish strangelets from MEMOs. Both forms of hypermatter would 
be extremely interesting to study and the discovery of one or the other would be worth 
every effort. Therefore, experiments should be devoted to the search for short-lived 
(anti-) hyperclusters. 

Large theoretical uncertainties remain, e.g. how the predicted yields depend on 
the model parameters. Both, theory and the current experimental results point towards 
a future search for sirange/eis/hyperclusters with rather short lifetimes. Experiments 
including a large TPC might be able to observe the decay short-lived hyperclusters. 
Indirect K + K~ correlation measurements might offer another possibility of detecting 
strangelets or hyperclusters | |249 . 



3. 7. Radiation of the plasma: direct photons and thermal dileptons 
Theoretical concepts 

The most prominent process for the creation of direct (thermal) photons in a QGP 
are qq — > jg (annihilation) and gq — > jq (Compton scattering). The production rate 
and the momentum distribution of the photons depend on the momentum distributions 
of quarks, anti-quarks and gluons in the plasma. Infrared singularities occuring in 
perturbation theory are softened by screening effects [ 346 , 347 , JJS, 349 1. If the 
plasma is in thermodynamic equilibrium, the photons may carry information on this 
thermodynamic state at the moment of their production |346| , [350j 351 , |352 . 



The main hadronic background processes to compete against are pion annihilation 
7T7T — > 7p and Compton scattering np — > 77T [ |346| , |353|| . The broad a\ resonance may act 
as an intermediate state in 7rp scattering and thus provide an important contribution 
354J , [353]1 via it's decay into 771-. In the vicinity of the critical temperature Tq a hadron 



gas was shown to "shine" as brightly (or even brighter than) a QGP [[346]. 

A finite baryochemical potential yields at constant energy density a reduced 
multiplicity of direct photons from a QGP ||355| , |356 . 



Hydrodynamical calculations can be used to compare purely hadronic scenarios 
with scenarios involving a first/second order phase transition to a QGP. They show 
a reduction in the temperature of the photon spectrum in the event of a first order 
phase-transition ||357| , |358| , |359|| . 
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The rapidity distribution of direct hard photons reflects the initial rapidity 
distribution of the produced mesons or directly the QGP ||360|| . It may thus provide 
insight into the (longitudinal) expansion of the photon source: If the hot thermal source 
is initially at rest and is accelerated by two longitudinal rarefaction waves propagating 
inwards with the velocity of sound, the photon rapidity distribution is strongly peaked 
around midrapidity. In contrast, a Bjorken-like boost-invariant expansion results in a 
more or less flat photon rapidity spectrum. 

If a very hot plasma is formed (e.g. at RHIC or LHC energies) a clear photon 
signal might be visible at transverse momenta in the range between 2 and 5 GeV/c 
[|361| , |362| , |363|1 - The lower p t range (1-2 GeV/c) is dominated by the mixed phase; 
separated contributions of the different phases are difficult to see due to transverse flow 
effects 1 13571| . These effects, however, can be important up to transverse momenta of 



5 GeV. Transverse flow effects also destroy the correlation between the slope and the 
temperature of the photon spectrum ||359| . 



Analogously to the formation of a real photon via a quark - anti-quark annihilation, 
a virtual photon may be created in the same fashion which subsequently decays into a 

pair (a dilepton). Also bremsstrahlung of quarks scattering off gluons can convert 
into dileptons. 

Dileptons can carry information on the thermodynamic state of the medium at 
the moment of production in the very same manner as the direct photons - since the 
dileptons interact only electromagnetically they can leave the hot and dense reaction 
zone basically undistorted, too. 

The main background contributions stem from pion annihilation, resonance decays 
364| |365| , |366| , |367| , |368|| (two pions can annihilate, forming either a virtual photon or a 



rho meson - both may then decay into a dilepton) and ir — p interactions | |369j , |370H at 
low dilepton masses and Drell-Yan processes [ 371 , 372] at high masses . Furthermore 
meson resonances such as the rho-, omega- or phi- meson may be produced directly 
or in the decay of strings and heavier resonances. As all of those vector mesons carry 
the same quantum numbers as the photon, they may decay directly into a dilepton. 
Resonances can also emit dileptons via Dalitz decays. The Drell-Yan process describes 
the annihilation of a quark of one hadron with an anti-quark (in proton proton collisions 
from the sea of q) of the other hadron, again resulting in a virtual photon which decays 
into a dilepton. The open charm contribution to the dilepton mass spectrum has been 
estimated to be negligible for low dilepton masses ||373|| at the CERN/SPS. At RHIC 
and LHC energies, however, charm contributions dominate the dilepton mass spectrum 



above 2 GeV p74 



Most original calculations on dileptons as signals of a QGP at CERN/SPS energies 



focused on masses below the rho meson mass ||375| , |376| , |377| , |378| , |379| , [I 
38 1| , |382] , |383|| . The current understanding of hadronic background contributions 
364| |366| , |367| , |368| | shows that most probably dileptons originating from a QGP are 



over-shined by hadrons, with the possible exception of masses around 1 to 1.5 GeV 



384J, [385|| where the rates from a plasma (at very high temperatures around 500 MeV) 



38 



may suffice to be visible. At higher masses, the yield of Drell-Yan processes from 
first nucleoli nucleon collisions most probably exceeds that of thermal dileptons from a 
QGP. Finite baryochemical potential will, at a given energy density, reduce the number 
of dileptons emitted from a QGP ||386|, |387|, |388| , due to the dropping temperature in 



that system. 

The dependence of the yield of high mass dileptons on the thermalization time is 
still a point of open debate [j389| , |390|1 . The parton cascade |391fl and other models of 
the early equilibration phase [|390| , |392| predict an excess of dileptons originating from 
an equilibrating QGP over the Drell-Yan background in the mass range between 5 and 
10 GeV. Then the early thermal evolution of the deconfined phase could be traced in 
an almost model independent fashion ||361|| . 

The secondary dilepton production via quark-antiquark annihilation has also 
been studied on the basis of a hadronic transport code (UrQMD [|110|| ). Here, one 
obtains a realistic collision spectrum of secondary hadrons for SPS energies. Using 
parton distribution functions and evaluating the contributions of all individual hadronic 
collisions one finds that meson-baryon interactions enhance the mass spectrum at mid- 
rapidity below masses of 3 GeV considerably ||372|| . Preresonance interactions are 



estimated to enhance this secondary yield by up to a factor of 5. 




Figure 19. Transverse momentum spectrum of directly produced photons in Pb+Pb 
collisions at 160 GeV/nucleon calculated with the UrQMD model. The resulting 
spectrum is compared with different hydrodynamical calculations. In all models the 
processes 7rry i— > 7r7, up t— > ttj and 7T7T t— > fry are considered as photon sources. The 
figure has been taken from |40£|. 
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Experimental status 

Experiments to measure direct photons are carried out at the CERN/SPS by the 
WA80/WA98 and the CERES/NA45 collaborations. Whereas the final analysis of the 
WA80 collaboration for S+Au indicates a 5% photon signal over background (with 
a 0.8% statistical and a 5.8% systematical error) |93], [394], [395| the CERES/NA45 
collaboration did not report any direct photons, but sets an upper limit of 7% for the 
integrated excess an unconventional photon source might have in central S+Au collisions 
396j |397| , PI . The WA80 collaboration has reported upper limits for each measured kt 



bin which yields important information for constraining the initial temperature of the 
reaction zone. Within the reported systematic errors the results of the WA80 and the 
CERES/NA45 collaborations are compatible with each other ||399| . 



Dileptons can be measured at CERN in form of dimuons by the HELIOS3, NA38 
and NA50 |400| , 401) , 402 , |403|| collaborations and in form of electron pairs by the CERES 
collaboration [ 397|| . Dimuons exhibit an excess in AA collisions in the mass range 
0.2 < M < 2.5 GeV/c 2 up to the as compared to pp and pA collisions. For 

dielectrons an excess is observed in the low-mass region 0.2 < M < 1.5 GeV/c 2 , again 



relative to pp and pA collisions (c.f. figure pip . 



Discussion 

The Pb+Pb analysis on direct photons of the WA98 and NA45/CERES collaborations 
is in progress. Hydrodynamical calculations are only compatible with the S+Au data 
of WA80 if a phase transition with its cooling is taken into account [|404| , [358| , |359| , [405 



406], \i07\\ or if higher mass meson and baryon multiplets are included for the hadronic 
EoS. 

Microscopic hadronic transport models, however, are not constrained by the 
assumption of thermal equilibrium, in particular in the initial stage, and yield results 
compatible with hydrodynamical calculations without invoking a phase transition 
scenario ||408|| , as can be seen in figure [19]. They shows that preequilibrium contributions 
dominate the photon spectrum at transverse momenta above ~ 1.5 GeV. The 
hydrodynamics prediction of a strong correlation between the temperature and radial 
expansion velocities on the one hand and the slope of the transverse momentum 
distribution on the other hand thus is not recovered in a microscopic transport model 



|40J 



Apart from these ambiguities in the interpretation of the data, the main problem 
with regard to the direct photon signal is the extremely small cross section in a difficult 
experimental situation, since photons from hadronic decays generate a huge hadronic 
background. The strong and dedicated effort to improve the measurements will be 
continued, also at the more promising collider energies. 

Both, the dielectron as well as the dimuon data seem to be compatible with a 
hydrodynamic approach assuming the creation of a thermalized QGP [[40911 • Hadronic 



transport calculations are not able to fully reproduce the observed excess ||410| , [411] , [412 
However, at least part of the observed enhancement of lepton pairs at intermediate and 
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low masses might be either caused by the previously neglected source of secondary Drell- 
Yan processes ||314|| or by contributions of heavy mesons, such as the a x [ |413|| . A detailed 
discussion of the dilepton data and its theoretical implications will follow in conjunction 
with the discussion on chiral symmetry restoration in section 



3.8. Restoration of chiral symmetry: vector mesons in dense matter 
Theoretical concepts 

The dilepton signal due to the decay of vector mesons, in particular from the rho 
meson, is of great interest: In conjunction with the chiral symmetry restoration 
T7|, ^18|| the rho, omega and phi mesons (and heavier meson resonances, 
are expected to change their spectral function in the hot, high 



e.g. the ax, ai 

baryon density medium: the breaking of the chiral SU(3)l x SU(3)ji symmetry (an 
approximative symmetry of QCD) results in quark condensates (qq) in the QCD vacuum 
and a "Goldstone" boson, i.e. the pion. The dependence of (qq) on the temperature T 
has been studied in the framework of lattice QCD [ [419| | and chiral perturbation theory 
420 , |421|| . Up to (0.7 - 0.8) Tq (qq) remains nearly constant and then its absolute value 
decreases rapidly (see figure |2^) . The behavior of the quark condensate at finite baryon 
densities is described in a model independent fashion by the Feynman-Hellman-theorem 
§422| , [4231 • A model calculation of the dependence of (qq) on both, the baryon density 
pj po and temperature T, can be seen in figure |(] - the drop of (qq) with p and T is quite 



analogous to the temperature and density dependence of the nucleon effective mass in 
the a — uj model as noted in |[424| . 



|<qq> P ,T 




T [MeV] 



Figure 20. quark condensate ~(qq) as a function of temperature T and baryon 



density p/po- The figure has been adapted from [414]. 



The reduction of the absolute value of the quark condensate (qq) in a hot and 
dense hadronic environment might reflect itself in reduced masses of vector mesons 
4251 B |27[ m, |2|, |43(J [431], [432], f|33f. However, a lowering of the mass of the p- 
meson - most commonly referred to as "Brown-Rho scaling" ||431|| - is not synonymous 



with the restoration of chiral symmetry: It has been shown by employing current 
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algebra as well as PCAC that to leading order in temperature, T 2 , the mass of the 
p- meson remains nearly constant as a function of temperature ||434|| , whereas the chiral 



condensate is reduced 11420 



Restoration of chiral may manifest itself in different forms | 417|| : the masses of the 
p- and the a\ meson may merge, their spectral functions could mix - resulting in peaks 
of similar strength at both masses (and causing a net-reduction at the p-peak) - or both 
spectral functions could be smeared out over the entire mass range. 

Dileptons from the in-medium decay of such vector mesons with modified masses 
and spectral functions would point towards the restoration of chiral symmetry at a phase 
transition. 



Experimental status 

At the BEVALAC the DLS collaboration has measured dielectron pairs in proton 



induced reactions as well as in d+Ca, He+Ca, C+C and Ca+Ca reactions [435, 436, 437 



Their latest results [437] for pair masses M < 0.35 GeV/c in the Ca+Ca system 
show a larger cross section than their previous measurements ||436|| and current model 
calculations |[438| , [439|| , suggesting large contributions from tc° and r\ Dalitz decays. The 
cross section da/dM scales with Ap ■ At up to pair masses of M = 0.5 GeV/c. For 
larger masses the Ca+Ca to C+C cross section ratio is significantly larger than the ratio 
of A P ■ A T values. 

Unfortunately, there are no experiments capable of measuring dileptons in A+A 
collisions at ~ 10 GeV/nucleon. However, in central sulfur - gold collisions at the 
CERN/SPS an enhancement has been measured in the invariant mass spectrum of 
muon pairs relative to the normalized proton proton and proton nucleus data at 200 
GeV/nucleon taken by the HELIOS/3 and NA38 collaborations: While the pp and 
pA data seem well described by measured sources such as Drell-Yan, open charm and 
hadronic decays, there is an excess in AA observed in the mass range 0.2 < M < 2.5 
GeV/c 2 (for the J/\l/, a suppression of the peak relative to this background is observed, 
c.f. section |3.9|) \^00[ |401| , f402j| . Similarly, for dielectron pairs in S+Au an excess has 
been observed by the CERES collaboration ||397|| in the low-mass region 0.2 < M < 1.5 
GeV/c 2 , again relative to pp and pA collisions (c.f. figure |21]). As in the case 



of the dimuon excess, the pp and pA data can be well understood taking known 
hadronic sources into account. Data for Pb+Au confirm this low mass dielectron excess 

EM MU rai. 



Discussion 

When searching for chiral symmetry restoration, thermal dileptons (see section |3.7| ) 
would serve as background. Due to the dominance of hadronic decays in the mass range 
up to 1.5 GeV, however, they do not pose a serious problem for the measurement of 
vector meson properties. 

As already stated in section [T7j, both, the dielectron as well as the dimuon data 
at the SPS seem explainable in a hydrodynamic approach assuming the creation of a 
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Figure 21. Left: inclusive e + e~ mass spectra in 200 GeV/nucleon S+Au collisions as 
measured by the CERES collaboration [ 397 1 . The figures have been taken from [444]. 
The shaded area depicts hadronic contributions from resonance decays. The data are 
compared calculations based on a purely hadronic scenario |41 



447, 44? 



Right: UrQMD prediction and data [441] for Pb+Au at 160 GeV/nucleon 



thermalized QGP |409 



On the other hand, the same data can be reproduced in the framework of 
microscopic hadronic transport models incorporating mass shifts of vector mesons 



445, 446, 447, 411, 410 



However, even bare hadronic transport model calculations, without any mass shift 
included, miss only the data in the 400 to 600 MeV bins (by 2 to 3 standard deviations) 



410, 411 



Calculations evaluating in-medium spectral functions, due to the coupling of the 
p with nucleon resonances and particle-hole excitations, also achieve a satisfactory 
reproduction of the CERES data ||449| , without requiring a dropping p-mass. Recent 
data on the p t dependence of this phenomenon ||450|| indicate however that the 
enhancement is most pronounced at low p t . This precludes the dominance of p-wave 
baryon-resonance effects. Since the hadronic transport models did - up to now - neglect 
contributions, e.g. from current-current correlation functions [[45 If and from additional 
heavy meson resonances, it has yet to be determined whether partial restoration of chiral 
symmetry or Brown- Rho scaling is the only possible explanation of these interesting new 
data. 

On the experimental side the main problems lie with the low signal to background 
ratio of 1/10 and limited statistics (< 1000 lepton pairs in Pb+Au). In addition, the 
shape of the e + e— excess at 0.25 GeV < m ee < 0.6 GeV coincides with that of the 
background. Currently a strong effort by the CERES collaboration is underway to 
upgrade the experiment with a TPC. The resulting increase in statistics and resolution 
should help to verify or falsify some of the conflicting hypotheses on the origin of the 
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low mass dilepton enhancement. 

The latest DLS data, gives rise to speculations that the observed enhancement in the 
mass-range below the free p meson mass may be due to enhanced p meson production or 
a mass shift of the p in a dense hadronic medium - without the need for a deconfinement 
phase transition. However, all microscopic model calculations which have addressed this 
data sofar ||452|, [453j|, have not been able to provide a reasonable explanation within the 



frameworks which work so well for dilepton production at CERN. 

3.9. Quarkonia Suppression: Evidence for Deconfinement or Dynamical Ionization 



In 1986 Matsui and Satz proposed [454] that the suppression of heavy quarkonia-mesons 
could provide one of the signatures for deconfinement in QCD at high temperatures. 
The idea was based on an analogy with the well known Mott transition in condensed 
matter systems. At high densities, Debye screening in a quark-gluon plasma reduces 
the range of the attractive force between heavy quarks and antiquarks, and above some 
critical density screening prevents the formation of bound states. The larger bound 
states are expected to dissolve before the smaller ones as the temperature of the system 
increases. The ip' and Xc states are thus expected to become unbound just above T c , 
while the smaller ip state may only dissolve above ~ 1.2T C . Heavier bb states offer the 



same features as cc states, but require much shorter screening lengths to dissolve [[455 
The T(bb) state may dissolve only around 2.5 T c , while the larger excited T' could also 
dissolve near T c . 

In order to determine the magnitude of suppression, it is obvious that the initial 
production mechanism must be well understood. Charm quark- anti-quark pairs, cc, 



are produced in rare pQCD gluon fusion processes [ |456|| , (gg — > cc) with a cross section 



in pp reactions a cS ~ 10 pb at yfs = 20 GeV. In the rare events when a pair is formed, 
both the charm and anti-charm quantum numbers remain approximately conserved, 
and either the cc emerge from the reaction in hidden charm quarkonium bound states, 
J/ip(lSi : 3097), ip'(2S\ : 3686), Xc(l-Po,i,2 : 3500), . . ., or in continuum open charm 

states D(1869), D*(2010), Even though only about 1% of the cc pairs emerge 

in pp collisions as J/ip states, these vector hidden charm mesons are the easiest to 
measure because they are seen as sharp resonances on top of a broad continuum in 
the invariant mass spectrum of di-leptons. In contrast, open charm production is much 
harder to measure. Semi-leptonic open charm decay contributes to the continuum yield 
of dileptons mainly below the peak. 

Above M^, the Drell-Yan (DY) process (qq — > pp) begins to dominate the 
continuum yield. The great importance of DY is that the absence of strong final state 
interactions of the produced leptons makes it possible to compute the absolute DY 
cross section via pQCD. The nuclear number dependence of the cross section is then 
entirely determined by geometrical (Glauber) factors, T^^b) (neglecting small nuclear 
dependence of the structure functions). Here, Tab{^>) is the number of binary NN 
interactions per unit clXGcl clS db function of the impact parameter. The measured DY 
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yields thus provide an important constraint on the impact parameter range associated 
with specific centrality (E?) triggers used in the experiment. The comparison of the 
centrality dependence of the J /if) and DY cross section therefore provides a calibration 
tool to determine the magnitude of the suppression factor of charmonium in nuclear 
collisions. 

Great interest in this proposed signature arose when NA38 found the first evidence 
of suppression in light ion reactions. With the new preliminary Pb + Pb data of NA50 
(see next section) which have been reported to show "anomalous" suppression, it is 
especially important to review critically some of the competing dynamical effects that 
could forge this quark-gluon plasma signature. 

One of the main problems in the interpretation of the observed suppression as 
a signal for deconfinement is that non-equilibrium dynamical sources of charmonium 
suppression have also been clearly discovered in p + A reactions. The interaction of the 
pQCD produced cc pair with any QCD medium (confined or not) decreases significantly 
the probability of that pair to emerge in an asymptotic J/if)(lS) state. The observation 
of J /if) suppression in p + A is direct proof of this fact since the formation of an 
equilibrated quark-gluon plasma in such reactions is not expected. 

A phenomenological analysis of p+A data yields a dissociation cross section for both 
the J /if) and the if)' around 7.3 ± 0.6 mb [ [457] , |458| , |161| , [459|| . This finding is surprising 
since the transverse areas of those two mesons differ by more than a factor of two. This 
has led to the ccg color octet model interpretation of the (quantum) formation physics 
involved. In this model, it is assumed that the pair is formed in a small octet state 
accompanied by a soft gluon that can be easily stripped off as it propagates through a 
nucleus [ |462| , |463|| . This qualitatively accounts for the observed equal nuclear absorption 
cross sections since, as a result of time dilation, hadronization into the asymptotic J /if) 



or if)' states is delayed several fm/c at high energies [ 463| , |161 



A recent development is the calculation of the hard contributions to the 
charmonium- and bottonium-nucleon cross sections based on the QCD factorization 



theorem and the non-relativistic quarkonium model [464]. The calculated p+A cross 



section agrees well with the data. The non-perturbative contribution to the charmonium 
cross section dominates at CERN/SPS energies and becomes a correction at LHC The 
J /if) production in nucleus-nucleus collisions at the CERN/SPS can then be reasonably 
well described by hard QCD, if the larger absorption cross section of the x states that 
are predicted by QCD are taken into account. 

While the octet model used in conjunction with the Glauber geometrical model 
can account for the pA observation, the corrections to this eikonal picture of nuclear 
absorption extrapolated to AA reactions are, however, not yet under theoretical control. 
For example, since the ccg color octet - soft gluon state is not an eigenstate of non- 
perturbative QCD, its effective hadronic absorption cross section may vary within the 
relaxation time. Also possible pileup of matter and energy loss prior to the gluon fusion 
event are neglected. The importance of gaining better theoretical control of the nuclear 
absorption process is underscored by the fact that in central PbPb reactions about one 
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half of the NA50 observed factor of four suppression of J /if) is estimated to arise from 
such non-equilibrium (quantum) formation physics. 

The second major theoretical uncertainty in interpreting charmonium suppression 
is distinguishing dynamical "background" dissociation processes such as if) + p — > DD 
and if) + A — > A C D from transient partonic dynamical processes (non-thermal color field 
fluctuations) and from the sought after screening mechanism in the plasma phase of 
QCD matter. 

Purely hadronic dissociation scenarios have been suggested [j46 5| , [466j [467j [468j [469 



470|1 , which could with suitable parameters account for J /if) and if)' suppression without 



invoking the concept of deconfinement. These hadronic scenarios are referred to as 
comover models. Suppression in excess to that due to preformation nuclear absorption 
is ascribed in such models to interactions of the charmonium mesons with comoving 
mesons and baryons which are produced copiously in nuclear collisions. Unfortunately 
none of the required absorption cross sections are experimentally known and estimates 
are highly model dependent. A general criticism of comover suppression estimates is 
the use of overly simplified Glauber geometry and idealized boost invariant expansion 
dynamics for the produced particles. 

Studying the transverse momentum, p t , dependence of J /if) and if)' production may 
yield additional information concerning the nature of the <//■?/>- suppression mechanism 
|471| , [472|| . Two common scenarios have been considered: At sufficiently high p t final 



state interactions might disappear due to time dilation, while hadronic absorption effects 
should be similar for the J /if) and if)' . In a deconfinement scenario this idea suggests that 
J /if) is suppressed only for low transverse momenta |[473|. A second scenario assumes 



that J /if) and if)' acquire large transverse momenta through multiple elastic parton- 
parton collisions. Those multiple collisions, however, are most likely to occur in the 
high density QGP region. The consequence would be that high p t J / if) and if)' should 
be even more suppressed than those with low transverse momenta [[474j| . A purely 
hadronic scenario predicts an increase in the mean transverse momentum as a function 
of transverse energy for the heavy Pb+Pb system [|475|| . 

The above discussion emphasizes some of the uncertain theoretical elements in the 
interpretation of charmonium suppression as a signal for deconfinement. To isolate the 
final state interaction effects from initial state nuclear absorption, it has been proposed 
to combine all available data using a Glauber model geometric variable, "L" . Since this 
is so popular we review below how this variable is defined. 

The suppression of the J /if) production cross section in A + B collisions can be 
expressed as 

a(AB -> J/if)) = ABa(pp -> J/if)e- [a ^ poL{A ' B ' ET)] S co . (3) 

Here o~^n ~ 5 — 7 mb is an effective preformation nucleon absorption cross section, po 
is the ground state nuclear density, and L(A, B, Ej) is a measure of the mean nuclear 
thickness evaluated through the Glauber model 

ABe -Klh P0 L{A,B,E T )] Sco= J d 2 b p {ETjb) J dz J Ssdz' 
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x p A (s, z)ps(|b - s|, z')T co (h, s) 



g u ccN 



fdz"{6{z"-z)p A (b,z")+e(z'-z")p B (b-s,z")) 

In the absence of preformation absorption (<7^ = 0) and the absence of comover 
absorption (T co = 0), the above factor reduces to AB for untriggered data (P = 1). As 
shown in the next section this AB scaling (expected for any hard pQCD process without 
final state interactions) is observed to hold very well for DY pair production. 

In addition to the uncertainties associated with (§]), the assumed constancy of the 
density along the path in the nuclear overlap region (neglect of energy loss and density 
pile-up especially at moderate SPS energies) and the assumed space-time independence 
of the effective cross sections of the pre-hadronic cc configuration, other sources of 
theoretical uncertainties are evident: A major source of model dependence of L enters for 
triggered data through the transverse energy impact parameter distribution, P(Ex,b). 
The observed transverse energy, E?, depends on the details of the experimental geometry 
and materials and is particularly difficult to simulate in the multi-target system of 
NA50. Often this distribution is simply parameterized such that its integral over impact 
parameters reproduces the the observed global dN/dEx distribution. The effective 
length L(E T ) is computed with the above assumptions by setting T co = 0. 

The main advantage of defining L is that data from different AB systems and E? 
triggers can be combined in one plot. However, we emphasize that unlike Et, L is not a 
measured quantity and is model dependent. Therefore, interpretations of data plotted 
as a function of L should be viewed with great caution. In contrast, it is theoretically 
much better to plot production cross sections as a function of AB to combine minimum 
bias data. 

The comover absorption factor S co depends sensitively on the magnitude and 
time dependence of the local comoving density of partons or hadrons as well as 
their absorption cross section. In additional, feed-down processes associated with 
"charmonium chemistry" must be taken into account in that calculation. The final 
J /if) include contributions from radiative decay of higher mass charmonium states. In 
one estimate ||476|| , it was assumed that p(ip' — > if)) ~ 12% of the observed if) arise from 
radiative ip' decay, and p(x — > if)) ~ 30% from \ decay. Unlike the small if), the larger if)' 
and x states are expected to have significantly larger absorption cross sections a co {nLJ). 
Evidence for comover absorption of if)' has been claimed to be observed already in S + U. 
Neglecting, for illustration, nuclear absorption and the impact parameter variations , 
the comoving survival factor is thus of the generic form 



co (nLJ) J^f dTp co (r) 

Even if a co {ip{\Si)) = 0, S co — > 0.6 if the higher mass charmonium states are absorbed 



S co = £ p (nLJ)e~ ac ° (nLJ)j To dTp - {T} . (5) 

nLJ 



Often simple scaling assumptions are assumed ||476|| for the evolution of the 
comoving density, p co oc dE±/dy(l/TR 2 ). With reasonable variations of the unknown 
parameters above, excellent fits to the NA38 S + U were obtained. However, even with 
these parameters fixed, comparisons for different AB systems require further dynamical 
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assumptions. Especially important is the assumed A-dependence of the the comoving 
matter density. Linear Glauber models tend to fail to reproduce the larger suppression in 
PbPb. However, nonlinear connections between p co {tq) and Et have been demonstrated 



to be also compatible with the data ||160 



Experimental status 

Systematic measurements of J/if>, if)' and 2-5 GeV continuum processes (Drell-Yan 
processes, open charm decay, etc.) have been performed by the NA38 collaboration using 



proton, oxygen and sulfur beams at CERN R4771 , [478! , 4T9 , 480 1 . The first preliminary 



data on Pb + Pb — > J /if) was reported by the NA50 collaboration in 1996. The data 



analysis is still not complete but in a recent conference [481] all the data have been 



combined as a function of AB as shown in Fig. |22|. The first striking result in the top 
left frame is that the Drell-Yan yield of dileptons with mass > 2.9 GeV scale within 
20% as a^y oc AB over five orders of magnitude in that variable. The so called K ~ 2 
factor depends on the choice of the proton structure functions. In the top right frame 
the deviation of the J /if) production cross section from this simple scaling is obvious 
(data are depicted with open triangles). After rescaling data from different energies to 
200 AGeV, all but the PbPb data lie on what appears to be a universal curve, (AB) ' 92 . 
The fact that the p + A data and S + U data lie on the same curve suggests the common 
preformation physics interpretation discussed before. The minimum bias PbPb data 
point is about 25% below the curve extrapolated up to 208 2 . This is the so called 
anomalous suppression. 

In the lower left and right frames the centrality Ex dependence of the J /if) and 
DY yields are compared. (Note that recently the Et scale has been recalibrated to 
be ~ 10 — 20 smaller than shown here, but no final publication has appeared as of 
this writing). The low Et < 40 GeV are more peripheral collisions with geometries 
comparable to central S + U. For Et > 50 GeV corresponding to b < 8 fm central 
collisions, the shape of the DY and J/ip distributions appear different. It should be 
noted that the J/ty and Drell-Yan distributions are affected by systematical errors which 
do not affect the ratio between the two variables and therefore the direct comparison of 
theory to these distributions may be misleading. The ratio of J /if) to Drell-Yan yields 
for Pb+Pb as a function of E t is displayed in figure [23| (full circles). Within the Pb+Pb 
sample a discontinuity has been hinted at around E t = 50 GeV. New minimum bias 
data with highly improved statistics show a rather smooth increase with decreasing Et- 
However, the data situation remains unclear for both, the low and the high Et limit. 
The (p 2 ) does not increase anymore for E t > 100 GeV in Pb+Pb collisions. 



Finally in Fig. |23] the provocative "L" plot | |481| | is shown that suggests that a 



sudden increase in if) absorption occurs for effective nuclear depths above ~ 7 — 8 fm, 
while the larger if}' state is similarly absorbed in S + U and and Pb + Pb. 

Discussion 

The initial plasma interpretations of the J/ip suppression in light nuclear beam data in 
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rescaled to 200 GeV, within 0<y<l and lcos(9 cs l<0.5 
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Figure 22. Preliminary J/V> to Drell-Yan data from NA50 |48l|, |48j]. The top left 
frame shows that Drell Yan scale with A x B over five orders of magnitude. Note that 
the value of the K is actually 1.8 and not 2.3 which refers to GRV LV. The top right 
shows a suppression of J /ip reaching almost a factor of 4 in Pb + Pb collisions (open 
triangles). The additional 25% suppression of J/tp in Pb+Pb relative to extrapolation 
of p + A and 5 + U is referred to as anomalous. The full circles depict a UrQMD 



calculation [483] including comoving mesons whereas the open circles show the same 
calculation with only nuclear absorption. In the bottom left and right frames the 
number of ip and DY pairs observed as a function of the uncorrected NA50 transverse 
energy are shown. The kink in the bottom left frame at Et ^ 40 — 50 GeV shows the 
rapid onset of anomalous suppression in the impact parameter range estimated to be 
b - 8 fm. 



1987 have been reformulated as a result of extensive p+A data proving the importance of 
pre-formation absorption phenomena in confined QCD matter. In addition, the similar 
suppression of ip and ip' as a function of A have allowed the determination of the octet 
cc effective cross section. The light ion suppression pattern and the p + A data are 
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Figure 23. Ratio of J/ip to Drell-Yan as a function of the experimentally measured 
transverse energy for the system Pb+Pb at 160 GeV/nucleon. The data are taken 



from 481], the histogram is a UrQMD calculation [483|. Both, data and calculation, 
have been rescaled to 200 GeV/nucleon. 
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Figure 24. Preliminary data on J/-0/DY and ip' /DY plotted as a function of the 
Glauber model nuclear thickness parameter, L, discussed in text. This plot suggests a 
sudden extra suppression of J/tp when the effective nuclear depth increases to about 
8 fm. 



now commonly agreed to be consistent with hadronic (confined matter) dissociation 
scenarios [|457| , [466j |467| , |161|1 - As pointed out , however, the pre-resonance ccg — g state 
is not an eigenstate of QCD and the use of an effective cross section for it in Glauber 
models is not without ambiguity. 

The data by the NA50 collaboration with the Pb+Pb experiment has given rise 
to renewed speculation on an additional suppression and on the possible creation of a 
deconfined phase ||484| , \i85[ [486|| . Plotting the total J/ip over Drell-Yan cross section 
ratio as a function of collision system mass (Ap ■ A?) an additional suppression in the 
order of 20% is observed between S+U and Pb+Pb. However, when studying the cross 
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section as a function of E?, a discontinuity at E? > 50 GeV may only be inferred 
when comparing the experimental J / ip over Drell-Yan ratio vs. E? with a Glauber- type 
calculation (without QGP formation) employing the same comoi>er-density as in the 
S+U case. The same Glauber-type hadronic absorption models ||470|| are also capable of 
reproducing the lead data of NA50 if a higher comov er-density is employed for Pb+Pb 
than for S+U [IH7]. For low transverse energies current hadronic comov er-models have 
difficulties in fitting the data. Therefore the issue whether simple comover models 
describe the lead data is not yet fully settled. 

Hadronic transport model calculations incorporate the full collision dynamics 
and go far beyond the commonly used simplified version of the Glauber theory 
|488j g|| , M |483| . These transport model calculations are very sensitive to certain 
input parameters such as the formation time of the J/ip and the comovers. In a first 
round of early calculations the HSD transport model ||488|| was fully able reproduce the 
NA50 lead data while assuming a fixed formation time of 0.7-0.8 fm/c for both, J/ijj 
and comovers. However, the hadron density which causes the suppression in the model 



may be unreasonably high. The UrQMD model [110], however, uses for the comovers a 



variable formation time emerging from the Lund string fragmentation formalism (here 
the formation time depends on the hadron mass) and zero formation time for the J/i/j. 
The assumption of zero formation time is valid if the J/ip is considered as pre-resonance 
ccg—g state with a hadronic dissociation cross section of 7 mb. However, in this mode the 
UrQMD model did not reproduce the additional suppression of the Pb+Pb experiment 
|489|| . The question of formation time might be a central issue since in the color octet 
model the dissociation cross section is actually higher during the lifetime of the pre- 
resonance cc$ — g state |[491| , |484j| than after hadronization. Furthermore the amount 
of comover - charmonium interaction will crucially depend on the formation time (and 
thus the cross section) of the comovers and whether the comovers are allowed to interact 
within their formation time as "pre-formed" states (analogously to the charmonium 
"pre-resonance" states). 

Recently, a new HSD calculation studied the influence of strings on the J/tp 
suppression. The cc pairs were produced perturbatively and the influence of dissociation 
of those pairs by strings was taken into account by regarding strings as longitudinal 
geometrical objects with a specific transverse radius ||490|| . Good agreement with the 
data from NA38 and NA50 was found for a string radius of R s ~ 0.2 — 0.3 fm. A 
new UrQMD calculation employs a microscopic Glauber simulation for J/ip production 
and the full microscopic transport calculation for nuclear and comover dynamics as 



well as for rescattering ||483|| . The dissociation cross sections are calculated using the 
QCD factorization theorem [f464|| , feeding from ip' (5%) and x states (40%) is taken into 
account and the cc dissociation cross sections increase linearly with time during the 
formation time of the charmonium state. Using only nuclear dissociation yields, a far 
smaller suppression than seen in the data is achieved. However, if comovers are taken 
into account (a meson ~ 2/3<7 mtc / eori , the agreement between theory and data is impressive 
(see figure |23[). The strong dependence of these results on details, such as the treatment 
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of the formation time or the time dependent dissociation cross section, remain to be 
studied further. 

Quantum effects such as energy dependent formation and coherence lengths must 
be taken into account |[492|| before definite statements can be made with regard to the 
nature of the J/ip suppression. 

Whereas there exist techniques to calculate the J /if) nucleon cross section without 
using the vector dominance model ||493| , the size of the cc-comover interaction is still 
unclear. TJNAF and HERMES experiments may be able to address the question 
of interaction of spatially small configurations with the nuclear medium. Such a 
research will help to understand the interaction of the cc wave packet with comovers at 
comparatively low energies [ 494}] . 

Interpretations of the data based on plasma scenarios are also increasingly evolving 
away from the original Mott transition analog. For example, J /if) suppression due 
to large coherent color fields (strings/ropes) have been proposed ||490|| . Percolation 
of longitudinal strings of transverse radius ~ 0.25 fm have been proposed to explain 
the possible sudden drop of the J/if> yield at moderate impact parameters (low Ej) 
495] . It is clear that much work needs to be done theoretically to sort through the 



many competing dynamical models of charmonium absorption. The observed effects are 
among the most striking results found in heavy ion reactions and deserve the intense 
attention they now receive. 

On the experimental side, it would be especially important to map out carefully 
the A dependence in the intermediate mass range 30 < A < 200 to confirm if there is 
a discontinuity or rapid change in the mechanism for moderate nuclear depths as was 
hinted at by the Pb data. In addition the functional dependence of the J/^f yield as 
a function of E? needs to be clarified for the low (Et < 40 GeV) and high E? {Ej- > 
100 GeV) limit. The beam energy dependence would be valuable to know, given the 
very rapid variation in pp at present SPS energies and to study the suggested expansion 
of small wave packets |[464| . Finally, an independent experimental confirmation of the 
results is essential. Most likely only at RHIC will there be an independent experiment, 
PHENIX, capable of addressing this very important observable. While nominally 
running at y/s = 200 AGeV, RHIC will be able to approach SPS conditions down 
to ~ 30 AGeV. 

Of course one of the most important additional experimental checks would be the 
search for discontinuities in other observables such as the strangeness fraction and HBT 
radii for less central reactions with E? in the region suggested by NA50 where new 
physics may arise. Thusfar, the other experiments have concentrated on more central 
collisions and the observables have varied smoothly as a function of control parameters, 
A,E, ... 

In summary, a striking pattern of charmonium suppression has been discovered 
by NA38/NA50 at the SPS. The theoretical debate on its interpretation is far from 
settled, but great strides have been made in the past decade to refine concepts and 
models. The rapid change in the suppression could be the smoking gun of deconfinement, 
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but it is not likely to be due to simple Debye screening effect originally hoped for. 
Rather novel QCD dynamics in a non-equilibrium plasma (e.g. enhanced color field 
fluctuations in moderate frequency u ~ 0.5 — 1.0 GeV) may emerge as the final culprit. 
A goal of further theoretical work will not be to continue to try to rule out more 
"conventional" explanations, but to give positive proof of additional suppression by 
QCD-based calculations which actually predict the ^-dependence of the conjectured 
signature. Consistency tests and a detailed simultaneous analysis of all the other 
measured observables are needed, if at least the same standards as for the present 
calculations involving other signatures are to be held up. 



4. Summary and Outlook 



AGS and SPS Milestones 

In the last two years the heavy ion research at Brookhaven and CERN have succeeded 
to achieve the measurement of a wide spectrum of observables with truly heavy ion 
beams Au + Au and Pb + Pb. As these programs continue to measure with greater 
precision the beam energy, nuclear size, and centrality dependence of those observables, 
it is important to recognize the major milestones past thusfar in that work. Experiments 
have conclusively demonstrated the existence of strong nuclear A dependence of 



baryon stopping power |134j , |143| 
hadronic resonance production [ 149j 

collective (transverse, directed, and elliptic) flow of baryons and mesons both at 
AGS and SPS energies [23(J 



strangeness enhancement [[293j , |295| , |303|, |144 



meson inter ferometric source-radii [|253 , |256] , 260 , 261 , 262 , 263 , 264 



dilepton- enhancement below the p mass ||397| , [440 



anomalous J / ip and ip' suppression |[496| , [48 2j , |48 1|| 



These observables prove that high energy- and baryon density matter has 
been created in nuclear collisions. The global multiplicity and transverse energy 
measurements prove that substantially more entropy is produced in A + A than simple 
superposition of Axpp would imply. Multiple initial and final state interactions play a 
critical role in all observables. The high midrapidity baryon density and the observed 
collective radial and directed flow patterns constitute one of the strongest evidence for 
the existence of an extended period (At ~ 10 fm/c) of strong final state interactions. 
The enhanced ip' suppression in S + U relative to p + A again attests to this fact. 
The anomalous low mass dilepton enhancement shows that substantial in-medium 
modifications of multiple collision dynamics exists, possibly related to mass shifts or 
in-medium broadening of vector mesons. The systematics of the strangeness (and anti- 
strangeness) quantum number production shows that novel non-equilibrium production 
processes arise in these reactions. Finally, the centrality dependence of J /ip absorption 
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in Pb + Pb collisions hints towards the non-equilibrium nature of such reactions, but can 
also be seen - in the case of an anomalous centrality dependence - as indication that high 
frequency gluon modes may be excited in such reactions. Is this the sought after quark- 
gluon plasma that thusfar has only existed as a binary array of predictions inside teraflop 
computers? It is too early to tell. Theoretically there are still too many "scenarios" and 
idealizations to provide a satisfactory answer. And there are experimental gaps such as 
lack of intermediate mass A ~ 100 data and the limited number of beam energies studied 
thusfar. The field is at the doorstep of the next milestone: A + Aa.ty/s = 30 — 200 
AGeV due to begin at RHIC/BNL in 1999. 

At the AGS, where particle spectra already have transverse slopes T > Tc, the 
highest chances for the discovery of partonic degrees of freedom lie in the measurement of 
the collective flow excitation function and the search for novel strangelet configurations. 
The investigation of the physics of high baryon density (i.e. partial restoration of chiral 
symmetry via properties of vector mesons), for which the gold beam at the AGS would 
be ideal, are unfortunately not accessible due to the lack of experimental setups capable 
of measuring electro-magnetic probes in AA collisions. 

At the CERN/SPS new data on electro- magnetic probes, strange particle 
yields (most importantly multistrange (anti-)hyperons) and heavy quarkonia will be 
interesting to follow closely. Energy densities estimated from rapidity distributions and 
temperatures extracted from particle spectra indicate that initial conditions should be 
near or just above the domain of deconfinement and chiral symmetry restoration. With 
respect to HBT (meson interferometry) source radii, the matter is not yet adequately 
resolved (extensive and precise comparisons with hadronic and deconfinement model 
calculations have yet to be performed). Directed flow has been observed - a flow 
excitation function, filling the gap between 10 AGeV (AGS) and 160 AGeV (SPS), 
would be extremely interesting to look for the softest point of the QCD equation of 
state. An effort to perform experiments in this energy region at the SPS is underway. 
However, dedicated runs would be mandatory to really explore these intriguing effects 
in the excitation function. It is questionable, whether this key program will actually get 
support at the SPS. Also the excitation function of particle yield ratios (ti/p, d/p, K/ir...) 
and, in particular, multistrange (anti-)hyperon yields, may be a sensitive probe of 
changes in the physics of the EoS. Most intriguing, however, would be the search for 
novel, unexpected, forms of SU(3) matter, e.g. MEMOS, strangelets or even charmlets. 
Such exotic QCD mesonic and nuclear configurations would extend the present periodic 
table of elements into hitherto unexplored flavor dimensions. A strong experimental 
effort should continue in that direction. The current status concerning the phase- 
diagram of nuclear matter is depicted in figure 

We note that exotic non-equilibrium phenomena such as Disoriented Chiral 
Condensates (DCC) that could effect very small transverse momentum pion spectra 
and charge/neutral meson fluctuations should also be continued to look for. These are 
highly speculative and scenario dependent phenomena, but worth to be searched for 
along with the above exotic mesonic and nuclear states. 
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Figure 25. Phase diagram of nuclear matter, adapted from 49, 497|. The two dashed 
lines mark the location of the expected phase boundary at its level of uncertainty. The 
data points mark thermal freeze out parameters deduced from AGS and SPS data, 
taking flow into account. The arrows indicate how the freeze out conditions may be 
reached during the expansion of the fireball. The two gray points show thermal model 
fits to UrQMD predictions for heavy collision systems. 



Another intriguing possibility of observing the transition from nuclear matter to 



deconfined quark matter may lie in the timing structure of Pulsar spin-downs [ [498 
Pulsars contain a huge amount of angular momentum and rotational energy. The 
emission of electro-magnetic radiation and electron-positron pairs (over a time-span of 
millions of years) causes a reduction in the angular velocity and thus also a reduction in 
the centrifugal deformation of the Pulsar. Consequently the interior density of the Pulsar 
increases and may rise from below the critical density for a phase transition to above the 
critical value. In the resulting phase-transition a conversion from rather incompressible 
nuclear matter to highly compressible quark matter will take place (starting at the center 
of the Pulsar) which reduces the radius of the Pulsar and causes and anomalous decrease 
of the moment of inertia with decreasing angular velocity. This is superposed on the 
normal reduction of angular velocity due to radiation loss. In order to conserve angular 
momentum, the deceleration of the angular velocity may decrease or even change sign 
resulting in a Pulsar spin-up. The time-span in which this effect may be observable 
is estimated to be in the order of 10 5 years. Since the mean life of Pulsars is around 
10 7 years, 1% of the 700 currently known Pulsars may currently be undergoing this 
phase-transition. 

New facilities: RHIC and LHC 

RHIC will begin operation in 1999 with four detectors: two medium scale ones, 
BRAHMS and PHOBOS, as well as two large scale detectors, PHENIX and 
STAR. BRAHMS (Broad RAnge Hadron Magnetic Spectrometer) is a conventional 
spectrometer (adapted from the AGS program) with particle ID, covering the cm 
rapidity range to 4. PHOBOS is a two arm magnetic spectrometer which will be 
able to measure low p t charged hadrons and leptons at selected solid angles. 
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PHENIX is a large solenoid with a variety of multi-purpose detector arrays; its 
goal is the multiple detection of phase transition signatures via the measurement of 



hadrons, leptons and photons in the same central rapidity bin 1 499|| . Apart from the 



QGP signatures which are already discussed, the PHENIX experiment will also search 
for a disoriented chiral condensate (DCC) [] |77], |501| , |502 |. 



The main emphasis of the STAR (Solenoidal Tracker At Rhic) detector will be the 
correlation of many (predominantly hadronic) observables on an event-by-event basisQ 
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The great energy range and beam target range accessible with RHIC will allow a 
dedicated systematic search for the quark-gluon phase matter at energy density an order 
of magnitude above the transition domain. This occurs not only because the rapidity 
density of hadrons is expected to be 2-4 times larger than in central SPS collision, but 
also because pQCD dominated mini-jet initial conditions are finally reached at collider 
(-^/s ~ 200 AGeV) energy range. A whole class of new signatures involving hard pQCD 
probes (high px and jets) becomes available. 

At yet higher energies at LHC, QGP research efforts and planning are centered 
around the ALICE detector. Its design is similar to that of STAR. However, also dimuon 
arms (like in PHENIX) are planned. ALICE will be the only large scale heavy ion 
detector setup at LHC. At y/s ~ 5 ATeV even bottom quarkonia are copiously produced 
and transverse momenta twice as high (p± ~ 60 GeV/c) will be readily measurable to 
probe even deeper into the multiparticle dynamics in a QGP. 

For applications to nuclear collision observables, an extension of the QGP concept 
to non-equilibrium conditions is required. The popular use of simple fireball models 
may provide convenient parameterizations of large bodies of data, but they will never 
provide a convincing proof of new physics. Microscopic transport models are required 
that can address simultaneously all the observables and account for experimental 
acceptance and trigger configurations. Present work in parton cascade dynamics is 
based largely on analogy of transport phenomena in known abelian QED plasmas. A 
significant new feature of QCD plasmas is its ultrarelativistic nature and the dominance 
of (gluon) radiative transport. These greatly complicate the equations. The role of 
quantum coherence phenomena beyond classical transport theories has only recently 
been established within idealized models. Much further work will be required in this 
connection. The outstanding theoretical task will be the development of practical (vs. 
formal) tools to compute quantum non-equilibrium multiple collision dynamics in QCD. 
Recent work [505] along the lines of non-compact lattice formulations of gauge theories 

+ A first preliminary result was published on the DCC search at CERN/SPS by the WA98 collaboration 
J500| : At a 90% confidence level they rule out a DCC admixture of greater than 20% 
* Recently, the NA49 collaboration presented first measurements of event-by-event fluctuations at 
the CERN/SPS, which may have significant implications for the issue of thcrmalization and critical 
fluctuations near a phase transition ||503| . The measured disappearance of strong dynamical fluctuations 
in (p t ) suggest a high degree of equilibration or at least rescattering. The absence of non-Gaussian 
fluctuations furthermore may exclude the possibility that the system has been close to a phase- 
transition. 



56 



may provide one of the most promising avenues in that direction. As yet unrealized 
techniques utilizing supersymmetry and string theory should also be explored. 

Experiments and data on ultra-relativistic collisions are essential in order to 
motivate, guide, and constrain such theoretical developments. They provide the only 
terrestrial probes of non-perturbative aspects of QCD and its dynamical vacuum. The 
understanding of confinement and chiral symmetry remains one of the key questions at 
the end of this millennium. 
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